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Abstract. We present a new synthetic model to follow the evolution of a planetary nebula (PN) and its central 
star, starting from the onset of AGB phase up to the white dwarf cooling sequence. The model suitably combines 
various analytical prescriptions to account for different (but inter-related) aspects of planetary nebulae, such as: 
the dynamical evolution of the primary shell and surrounding ejecta, the photoionisation of H and He by the 
central star, the nebular emission of a few relevant optical lines (e.g. H/3; HellA4686; [O ill] A5007). Particular 
effort has been put into the analytical description of dynamical effects such as the three-winds interaction and 
the shell thickening due to ionisation (i.e. the thin-shell approximation is relaxed), that are nowadays considered 
important aspects of the PN evolution. Predictions of the synthetic model are tested by comparison with both 
findings of hydrodynamical calculations, and observations of Galactic PNe. The sensitiveness of the results to the 
models parameters (e.g. transition time, mass of the central star, H-/He-burning tracks, etc.) is also discussed. We 
briefly illustrate the systematic differences that are expected in the luminosities and lifetimes of PNe with either 
H- or He-burning central stars, which result in different "detection probabilities" across the H-R diagram, in both 
H/3 and [Olll] A5007 lines. Adopting reasonable values of the model parameters, we are able to reproduce, in a 
satisfactory way, many general properties of PNe, like the ionised mass-nebular radius relationship, the trends of 
a few main nebular line ratios, and the observed ranges of nebular shell thicknesses, electron densities, and ex- 
pansion velocities. The models naturally predict also the possible transitions from optically-thick to optically-thin 
configurations (and vice versa). In this context, our analysis indicates that the condition of optical thinness to the 
H continuum plays an important role in producing the observed "Zanstra discrepancy" between the temperatures 
determined from H or Hell lines, as well as it affects the mass- increasing part of the ionised mass-radius relation. 
These predictions are supported by observational indications by Mendez et al. (1992). Another interesting result 
is that the change of slope in the electron density-nebular radius relation at i?ion ~ 0.1 pc, pointed out by Phillips 
(1998), is also displayed by the models and may be interpreted as the result of the progressive convergence of the 
PNe to the condition of constant ionised mass. Finally we would like to remark that, thanks to its computational 
agility, our synthetic PN model is particularly suitable to population synthesis studies, and it represents the basic 
ground from which many future applications will be developed. 
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1. Introduction 

Planetary Nebulae (PN) are among the most beautiful 
and fascinating celestial objects. They appear when low- 
and intermediate-mass stars, after losing their external en- 
velopes during the AGB phase, cross the H-R diagram on 
their way to the white dwarf cooling sequence. During 
this crossing, the central nuclei become sufficiently hot to 
ionise the circumstellar material previously expelled. The 
resulting emission nebulae may turn bright enough (spe- 
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cially in the [O ill] lines) to emerge strikingly against the 
background of "normal" stars, even at relatively large dis- 
tances. 

The study of the spectra of PNe provide unvaluable 
information on the chemical enrichment of the insterstel- 
lar medium from low- and intermediate-mass stars, and 
on the kinematics of their parent populations. In general, 
they offer a unique tool to test stellar evolution theory 
(e.g. mixing processes, stellar atmospheres, evolutionary 
time-scales, etc.), as well the many physical processes in- 
volved (e.g. stellar winds, gas dynamics and kinematics, 
photoionisation, line formation, etc.) The reader can get 
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an overall picture by referring to the works of Pottasch 
(1984), Peimbert (1990), and Kwok (2000). 

Indeed, the observed properties of PNe are the result 
of the complex interplay between several factors, whose 
inter-relationship cannot be summarised in simple terms. 
For instance, the PN lifetime - i.e. the detectability period 
- crucially depends on both i) the evolutionary time-scale 
of the central star (the time it takes to become and re- 
main an efhcient emitter of ionising photons) and ii) the 
dynamical time-scale of the circumstellar ejecta (including 
the duration of the ejection event on the AGB tip, and the 
time taken by the nebula to expand away from the central 
star). 

It is clear that the detectability of a PN requires a suit- 
able tuning between the two time-scales. If the central star 
evolves too quickly the brightening of the ionised nebula 
would correspond to a short-lived event having a low de- 
tection probability. Conversely, if the central star evolves 
too slowly it would start to ionise the nebula when this 
has already dispersed to too large radii to be detected. It 
follows that only when these two time-scales are compara- 
ble, we can have a maximum probability of detecting the 
PN phase. As a matter of fact, neither the stellar evolu- 
tionary time-scale nor the dynamical time-scale are simple 
predictable functions of the involved parameters. 

Though all post- AGB stellar models indicate that, in 
general, the evolutionary speed on the H-R diagram de- 
creases at increasing final stellar mass, other factors play 
non-negligible roles (such as: past stellar history, phase 
of the pulse cycle marking the termination of the AGB, 
efficiency of the fast wind, chemical composition of the 
envelope; see Vassihadis & Wood 1994; Blocker 1995 for a 
detailed discussion). This complexity translates into quite 
large uncertainties affecting, in particular, the so-called 
transition time (that determines the onset of the photoion- 
isation phase), and the fading time (that determines the 
end of the PN lifetime). 

The dynamical time-scale of the nebula is also deter- 
mined by several factors including i) the properties (mass- 
loss rates and terminal velocities) of both the AGB (super- 
and/or slow) wind and the post- AGB fast wind, ii) their 
dynamical interaction producing gas shocks with related 
heating/cooling processes, and iii) the effect due to ion- 
isation which tends to expand and accelerate the nebula 
as a consequence of the increased thermal pressure. 

The importance of these dynamical factors has been 
pointed out by several past analytical works (e.g. Kahn 
1983, 1989; Volk & Kowk 1985; Kahn & Breitschwerdt 
1990; Breitschwerdt & Kahn 1990), and then broadly con- 
firmed and fully analysed by (magneto-)hydrodynamical 
simulations of PN evolution (e.g. Schmidt- Voigt & 
Koppen 1987ab; Marten & Schonberner 1991; Balick et 
al. 1992; Mellema & Frank 1995; Schonberner & Steffen 
2000, and references therein). 

If, on one hand, hydrodynamical simulations are able 
to solve the fine details of the nebular structure get- 
ting insight into the involved physical processes, on the 
other hand they are so demanding in terms of comput- 



ing time that only few cases of PN evolution have so far 
been adequately modelled. The hydrodynamical approach 
does not permit a systematic study of PN properties as 
a function of the fundamental parameters yet, i.e. stellar 
mass and metallicity, transition time, properties (mass- 
loss rates and terminal velocity) of stellar winds, etc. As a 
consequence, it becomes difficult also to predict the global 
properties of PNe populations in different galaxy environ- 
ments - e.g. how they should appear for varying histories 
of star formation and chemical enrichment. 

These aims may be achieved, taking a reasonable com- 
putational effort, with the aid of "synthetic" models of 
PNe, in which the main factors determining their evolu- 
tion are described by means of reasonable - and neces- 
sarily simple - approximations. Several examples of syn- 
thetic PNe models can be found in the literature, as in e.g. 
Kahn (1983, 1989), Kahn & West (1985), Volk & Kowk 
(1985), Stasihska (1989), Kahn & Breitschwerdt (1990), 
Breitschwerdt & Kahn (1990), Stasihska et al. (1998), 
Stanghellini & Renzini (2000). 

However, it is worth making a due distinction within 
the class of synthetic models. In fact, one can opt for 
either i) schemes based on rather crude assumptions 
(e.g. constant nebular mass and/or expansion velocity) 
which do not attempt to predict the dynamical properties 
and/or the ionisation structure of PNe (e.g. Stanghellini & 
Renzini 2000), or ii) more detailed analytical approaches 
which account for an overall, even if simplified, description 
of PN evolution, including the dynamical one (e.g. Volk & 
Kowk 1985; Kahn & Breitschwerdt 1990). 

As far as the dynamical aspect is concerned, a further 
distinction should be made. Most dynamical models for 
PNe (both synthetic and hydrodynamical) may be classi- 
fied as either Three-Wind Models if they consider the pos- 
sible interaction between the AGB slow wind, the AGB 
superwind, and the post-AGB fast wind (e.g. Schmidt- 
Voigt & Koppen 1987ab; Marten & Schonberner 1991), or 
Two-Wind Models if they lack of a distinguishable super- 
wind phase (e.g. Kwok et al. 1978; Volk & Kwok 1985), 
whose existence has actually been supported by observa- 
tional evidence (see, for instance, van Loon et al. 1998). 

The present study belongs to the class of synthetic 
Three-Wind models for PNe and accounts, as far as possi- 
ble, for many relevant aspects, including the dynamics of 
the nebula, the ionisation structure of H and He by adopt- 
ing stellar fluxes from non-LTE model atmospheres, and 
the emission of a few recombination and forbidden lines. 
Of course, we cannot avoid to adopt some simplifying ap- 
proximations and assumptions - like spherical symmetry, 
constant electron temperature. However, we will show that 
several tests (both theoretical and observational) support 
the overall performance of the model, so that it can be 
considered as a satisfactory tool for investigating the gen- 
eral properties of PNe. 

The outline of the paper is as follows. The basic ingre- 
dients of the model are fully described in Sect. ^. Then, 
after summarising the basic assumptions of the models, 
in Sect. ^ we give a brief overview of the model predic- 
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tions. We proceed in Sect. ^ with the exploration of the 
predicted properties of PNe as a function of the stellar 
mass, and of some free model parameters (e.g. transition 
time, electron temperature, H-/He-burning track). To this 
aim, we consider several observables to be reproduced, al- 
ways limiting to samples of Galactic PNe. In this way, as 
a first step, we limit the possible effects due to different 
metallicities, possibly present in PNe populations of dif- 
ferent galaxies. Concluding remarks and future aims are 
expressed in Sect. ||. 

2. A new synthetic model of PN evolution 

2.1. The general scheme 

Before proceeding with the detailed description of the 
model, it is convenient to sketch the general scheme of 
the several different processes here considered. 

— We start with detailed synthetic models of AGB stars 
of varying initial mass and metallicity. According to 
the adopted prescription, mass loss gets more and more 
efficient as the star climbs the AGB, so that most of the 
stellar envelope is ejected during the last "superwind" 
stages. We attribute proper expulsion velocities to the 
AGB ejecta, and follow their expansion during all the 
subsequent evolution. 

— The central star, peeled off of its envelope, is assumed 
to evolve along the post-AGB track determined by its 
mass and initial metallicity. The so-called "transition 
time" , corresponding to the very initial post-AGB evo- 
lution, is left as a free parameter. At each time step, 
we compute the "fast wind" (momentum and kinetic 
energy) emitted by the central star, as well as the flux 
of ionising photons. 

— The fast wind shapes a shell nebula through its inter- 
action with the slower stellar wind emitted during the 
AGB evolution. We compute the resulting dynamical 
evolution by means of the so-called "interacting- winds 
model" . As the shell expands, it continuously accretes 
matter from the circumstellar layers, ejected at earlier 
ages and with lower velocities during the AGB phase. 

— As soon as the flux of ionising photons becomes large 
enough (due to heating of the central star), an ioni- 
sation front proceeds through the nebula. We follow 
its position as a function of time, as well as the to- 
tal emitted fluxes in a few relevant recombination and 
forbidden lines. 

— The additional pressure provided by ionisation pro- 
duces a thickening of the shell, that we follow in an 
approximative way. 

— Eventually, the nebula expands to very low densities, 
and the central star fades along its white-dwarf cooling 
sequence, causing both dispersion and dimming of the 
emission nebula. We follow the evolution of all nebular 
quantities up to sufficiently large times. 

The subsections below will illustrate in detail all the pre- 
scriptions we have adopted. 



2.2. The AGB phase 

The AGB phase is described by detailed synthetic models, 
covering the evolution of low- and intermediate-mass stars 
from the first thermal pulse up to the complete ejection of 
the envelope (details can be found in Marigo 1998, 2001; 
Marigo et al. 1996, 1998, 1999). 

From these models we derive the mass of the central 
star as a function of the initial mass and metallicity of 
the progenitor, the mean radial gradient in density and 
chemical composition of the ejected material as a function 
of time (assuming spherical symmetry and stationary wind 
flows) . A simple attempt to account for possible dynamical 
interaction of AGB shells ejected at different epochs with 



different velocities is presented in Sect. 2.4.3 



Mass-loss on the AGB is included according to the 
semi-empirical formalism developed by Vassiliadis & 
Wood (1993; hereinafter also VW93), which combines ob- 
servations of variable AGB stars (Mira and OH/IR stars) 
with standard predictions of pulsation theory. 

In brief, mass-loss during the AGB consists of two 
phases. Initially, for P < 500 days, the mass-loss rate ex- 
ponentially increases with the pulsation period P 



logAf = -11.4 + 0.0123 P 



(1) 



until when, for longer periods, it attains and maintains 
nearly constant values, typical of the super-wind regime 



M = 6.07023 X 10-^/3 



cVagb 



(2) 



where f3 = 1.13 (Z/0.008) following Bressan et al. (1998). 
Equation (^) is derived simply assuming that all the mo- 
mentum carried by the stellar radiation fleld is transferred 
to the outer gas, and allowing the possibility for multiple 
scattering of photons (/? > 1; typically /3 ~ 2 for solar 
metallicity) . 

In Eqs. (|^) and (^, M is given in units of Mq yr^^, the 
stellar luminosity L is expressed in Lq, the pulsation pe- 
riod P in days, the light speed c and the terminal velocity 
of the stellar wind Vagb are given in km s^^. 

The pulsation period P is derived from the period- 
mass-radius relation (equation 4 in VW93), with the as- 
sumption that variable AGB stars are pulsating in the 
fundamental mode (see Wood et al. 1999): 



logP = -2.07+ 1.94 logP- 0.9 logM 



(3) 



where the period P is given in days; the stellar radius 
R and mass M are expressed in solar units. The wind 
expansion velocity Vagb (in km s"""^) is calculated as a 
function of the pulsation period P (in days): 



Vagb = -13.5 + 0.056 P. 



(4) 



In VW93 the allowed variability range was constrained 
according to observations of Mira and OH/IR stars. 
Specifically, lower and upper limits of V'agb ~ 3km s~^ 
and V^QB = 15km s^^, respectively, were adopted. 

At this point it should be remarked that, for our pur- 
poses, V^AGB is ^ quite crucial quantity, since it essentially 
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determines the characteristic velocity of the superwind 
ejecta, from which most PNe are thought to originate. 
Therefore, setting V^cb practically equivalent to spec- 
ify the expansion velocity of the main ejection event at the 
end of the AGB phase, which is clearly of great relevance 
to the subsequent dynamical evolution of the PN. A few 
alternative prescriptions for V^^g are presented below. 

2.2.1. The maximum velocity 

Several factors affect the terminal velocity of the AGB 
wind. Observations of OH/IR and Mira variable stars in- 
dicate that Vagb may range over a relatively large range, 
i.e. from 8 to 30 km s~^, with a typical mean value 
of ~ 14 — 15 km s^^ for galactic variables (e.g. Gussie & 
Taylor 1994; Groenewegen & de Jong 1998; Groenewegen 
et al. 1998, 1999; van Loon 2000), and possibly lower 
for variables belonging to more metal poor galaxies (e.g. 
Vagb (LMC)/ Vagb (Galaxy) ~ 0.5 - 0.6 for OH/IR stars, 
as pointed out by Wood et al. 1992). In the framework of 
radiation-driven dusty winds, this fact would be the re- 
sult of lower dust-to-gas ratio, 5, at decreasing metallicity 
(e.g. Vagb oc S°'^ according to the theoretical expecta- 
tions of Habing et al. 1994). Moreover, both observations 
(e.g. Loup et al. 1993) and model predictions (e.g. Netzer 
& Elitzur 1993) converge to the result that Vagb is larger 
for carbon-rich stars than for oxygen-rich stars, owing to 
the different absorption properties of chemically different 
dust grains. 

In addition to a dependence on chemical composition, 
calculations of radiative transfer in dust-driven winds (e.g. 
Netzer & Elitzur 1993; Habing et al. 1994; Kriiger et al. 
1994; Ivezic & Elitzur 1995) have pointed out that the 
terminal velocity is affected by several other factors, i.e. 
mass-loss rate, radius of dust grain, and stellar luminosity. 

A theoretical finding, which is relevant to the present 
analysis, is that at sufficiently high mass-loss rates (M 
1O~^M0), the terminal velocity is expected to be practi- 
cally independent from M (Habing et al. 1994). 

On the basis of all these aforementioned elements 
(both theoretical and observational), let us now consider 
some prescriptions we may reasonably adopt for V^Qg, 
namely: 



1. 



V^GB = constant 



say 15 km s~ , regardless of any 
other parameter. This is the simplest possible assump- 
tion, and is made for instance, by VW93 and Marigo 
et al (1996, 1998). 

V^(^^ = f{Z) = 6.5 (Z/0.008) + 0.00226 P , (5) 
which includes a metallicity dependence, according to 
an empirical calibration based on OH/IR variables of 
the Galaxy and LMC (Bressan et al. 1998). Such a 
prescription has been employed in the AGB synthetic 
calculations adopted here (Marigo et al. 1999; Marigo 
2001). For typical values of P > 500 - 2000 days, we 
get V^^g - 8 - 11 km s'^ for Z = 0.008, 



and V^gI 



30 



20 - 



10 



4.2 
log (L/LJ 



4.4 



4.6 



Fig. 1. Expected dependence of Vagb on stellar lumi- 
nosity, according to Kriiger et al. 1994 (squares). The 
parabohc fitting relation (Eq. ^) is shown (sohd hue), to- 
gether with that obtain ed sca ling it by a factor = 0.588 
(dotted hue). See Sect. ^.2.1 for explanation 



metallicity, V^'gb almost constant, due to the weak 
dependence on the pulsation period. 

=fiL)=A (log L)^ + B\ogL + C, (6) 
which includes a luminosity dependence, according to 
the results of two-fluid-model calculations for dust- 
driven winds (Kriiger et al. 1994). The fitting coef- 
ficients (to reproduce their figure 12, see also Fig. ^ 
are A = -16.7282, B = 164.1858, and C = -366.0314. 
Of course, these numbers only apply to a specific set 
of calculations performed by Kriiger et al. (1994). 
Therefore, in order to specify V^gb model, it 

may be reasonable to keep the relative trend with the 
stellar luminosity as predicted by Kriiger et al., and 
scale it by a proper multiplicative f acto r, /„, calibrated 
on some suitable constraint (Sect. 4.1). 



17 - 20 km s^i for Z = 0.02. We note that, at fixed 



2.3. The post-AGB evolution of the central star 

The post-AGB tracks are taken from Vassiliadis & Wood 
(1994) data set, that spans a large range of core masses 
0.55 Mq - 0.95 Mq) and metallicities {Z =0.016, 0.008, 
0.004, and 0.001). They also include both H- and He- 
burning central stars, i.e. stars that have left the AGB at 
different phases of their thermal pulse cycles, while they 
were efficiently burning their nuclear fuel either in the 
H- or in the He-burning shell. These two situations lead 
to very different lifetimes during the post-AGB, including 
the possible presence of a late flash and "final loop" for 
He-burners. All evolutionary tracks start at log Toff — 4 
which, in the case of He-burning stars, marks the point in 
the blue-ward part of the last loop in the H-R diagram. 
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Fig. 2. Evolutionary tracks in the H-R diagram from the 
first thermal pulse on the AGB to the cooling sequences 
of White Dwarfs. The transition from the tip of synthetic 
AGB models to the starting point of the post-AGB H- 
burning tracks of the same core mass is shown (dotted 
line) 

For all the details about the models the reader should refer 
to the Vassiliadis & Wood (1994) paper. 

At fixed metallicity, the match between the adopted 
AGB synthetic models and the post-AGB tracks is per- 
formed by interpolating the latter data as a function of 
the stellar mass left at the end of the AGB. An example 
is illustrated in Fig. ||. 

Minor differences between the stellar luminosity at the 
tip of the AGB and that of the horizontal part of the 
post-AGB track in the H-R diagram reflect small differ- 
ences between the core mass-luminosity relation derived 
by Vassiliadis & Wood (1993, 1994), and that adopted in 
the synthetic AGB models here in use (Wagenhuber & 
Groenewegen 1998). 

As we can see, the early post-AGB stages (i.e. after the 
star has left the AGB up to log Teff = 4) are not covered by 
the tracks (dotted line in Fig. |2|). They correspond to the 
transition from the end of the AGB to the starting of the 
ionisation phase, i.e. the period in which the central star 
does not emit enough UV photons to ionise the H atoms 
in the surrounding nebula. Large uncertainties affect the 
duration of such phase - the so-called transition time, tu 
(see Renzini 1989; Schonberner 1990) - essentially due to 
our poor knowledge of the properties of the stellar wind, 
as it evolves from the super-wind regime to the fast-wind 
regime (see Blocker 1995 for an analysis of this point). 

Moreover, even for a given prescription for mass-loss, 
Vassiliadis & Wood (1994) find no evident correlation be- 
tween ttr and the mass of the central star, but rather a 



dependence on the phase of the pulse cycle at which the 
star leaves the AGB (Vassihadis & Wood 1994). For all 
these reasons, we choose to adopt itr as a free parameter 
in our model, comprised in the range 500 yr - 10 000 yr. 

2.3.1. The fast wind 

Mass-loss during the post-AGB phase of the central star is 
a crucial process as it affects the evolutionary time-scales 
and hence observable properties of PNe. 

In this model, we describe the fast wind adopting the 
prescriptions from Vassiliadis & Wood (1994). This for- 
malism combines the results of the radiation-pressure- 
driven- wind theory for hot stars (Pauldrach et al. 1988), 
with empirical correlations observed for both AGB stars 
and central stars of planetary nebulae. In this way, a con- 
tinuous change of the wind parameters is allowed, as the 
evolution proceeds from the super-wind regime, typical of 
the last stages of the AGB phase, to the accelerating fast 
wind regime during the post-AGB phase. 

Let us summarise the basic prescriptions. The terminal 
wind velocity, v, of the fast wind is calculated with 



log = -2.0 + 0.52 log Teff, 

^CSC 



where the escape velocity is 



(i-r) 



2GM 
R 



1/2 



(7) 



(8) 



and r is the ratio of the stellar to the Eddington luminos- 
ity. This latter is given by 



^Edd — 



0.2{l + X) 
inc GM 



(9) 



where the other symbols have the usual meaning. 
Equation (0) is a linear fit to observed data, including 
central stars of PNe, B stars, A-F supergiants, and a fidu- 
cial point for OH/IR stars. 

Let us denote by Mii,„ the radiation-pressure-driven 
mass- loss limit, corresponding to the case in which all the 
momentum carried by the radiation field is transferred to 
the wind gas: 



lim 



L 

cv 



(10) 



Then, the mass-loss rate, M, in the PN phase is calcu- 
lated assuming a linear relation between log M/Mn^ and 
log Teff: 



log = -3.92 + 0.67 log Teff 

Mini 



(11) 



as derived by Vassiliadis & Wood (1994) following the the- 
oretical results by Pauldrach et al. (1988). 
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Fig. 3. Schematic sketch of the system composed by the central star, shocked PN shell, and undisturbed AGB wind 
2.4. The dynamical evolution of the shocked nebula Table 1. Adopted notation for relevant quantities 



The dynamical evolution of planetary nebulae is followed 
in the framework of the interacting-winds model, first de- 
veloped by Kwok et al. (1978), Kwok (1983), Kahn (1983), 
and substantially improved by Volk & Kwok (1985, here- 
inafter also VK85). Major improvements to the VK85 
scheme are included, by combining complementary pre- 
scriptions from other studies (i.e. Kahn 1983; Kahn & 
Breitschwerdt 1990; Breitschwerdt & Kahn 1990), as de- 
scribed in the next sections. 

For the sake of clarity, only the basic concepts and 
definitions are briefly recalled in the following (see also 
Table 1). For more details the reader should refer to the 
aforementioned papers. 

The dynamical evolution of PNe is considered the re- 
sult of the interaction between the fast wind emitted 
by the central star during the post-AGB evolution, and 
the slow circumstellar wind previously ejected during the 
AGB phase. When the fast-moving gas catches up with 
the slower one, a dense shell forms and two shock waves 
are generated. The inner shock is located at the point 
where the fast wind encounters the slow- moving gas. The 
outer shock is at the outer boundary of the shell, beyond 
which the undisturbed AGB wind is expanding outward. 
Depending on the efficiency of radiative cooling, the hot 
shocked gas from the fast wind may fill (or not) the bubble 
extending from the inner shock up to the so-called contact 
discontinuity, i.e. the inner rim of the shell. No heat or gas 
is assumed to pass across it. The shocked shell may be par- 
tially or totally ionised by the UV photons emitted by the 
central star. 

According to the dominant mechanism causing the ex- 
pansion of the shell, we can distinguish two dynamical 
phases, namely: the initial momentum-driven phase and 
the subsequent energy-driven phase. 



Magb mass-loss rate during the AGB phase 

m mass-loss rate during the post-AGB phase 

Ms mass of the shell 

_Ri radius of the inner shock 

Rs radius of the contact discontinuity (thin-shell) 

Rq radius of the contact discontinuity (thick-shell) 

Rion radius of the ionisation front (thick-shell) 

7?out radius of the outer shock 

Vagb terminal velocity of the AGB wind 

V terminal velocity of the fast wind 

Vi velocity of the inner shock front 

Vs expansion velocity of the shell (thin-shell) 

Vg* velocity related to -Rs* (thick-shell) 

Vion velocity of the ionisation front (thick-shell) 

P pressure of the hot shocked gas 



Hereinafter we convene to set the evolutionary time 
equal to zero {t = 0) at the end of the AGB phase, i.e. at 
the onset of the momentum-driven phase. 

2.4.1. The momentum-driven phase 

During the initial stages after the end of the AGB phase, 
the rate of radiative cooling is larger than the rate of ki- 
netic input from the fast wind (mw^/2), so that the inner 
shock can be considered isothermal. The shocked matter 
can be compressed effectively and the inner front is pushed 
quite close to the contact discontinuity (i?i w Rs), as it 
turns out from hydrodynamical calculations (e.g. Marten 
& Schonberner 1991). Under these conditions, the driving 
force which makes the shocked shell move outward is the 
transfer of momentum from the fast wind. In our model, 
we account for this phase by adopting the prescriptions 
presented in Kahn & Breitschwerdt (1990, hereinafter also 
KB90). 
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is set to zero at the end of the AGB 
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From the conditions of mass and momentum conser- 
vation applied to the system composed by the fast wind 
and the swept shell, one gets simple relations expressing 
the time evolution of the shell mass 



where g = 4 x 10^ cm° g ^ s ^ (Kahn 1976) and all other 



dMs 
"dT 



MAGB(i) 



V^s(t) 



1 



(12) 



V^{t) ^ 1 + vM^) ^agbW, 



(13) 



yAGB(t) 

and the radius of the contact discontinuity 
di, 

where a is defined according to: 

m[t)v{t) = a(i)MAGB(i)^AGB(0 (14) 

For typical values of the winds parameters involved (e.g. 
Magb ^ 10-"^ - IO-^Mq yr-i, Vagb ~ 10 - 20 km s-\ 
rh ~ 10""^ - IO^^Mq yr-\ w ~ 100 - 1000 km s^^) the 
quantity a is found to be less than unity (~ 0.01 — 0.1). 
This means that the rate of momentum injection to the 
gas from the accelerating fast wind is lower than that due 
to the AGB superwind, i.e. mv < A^agb^gb- 

Under the condition of a steady rate of momentum 
injection, i.e. if a = const., the shell's motion would be 
uniform, as assumed by KB90. Here we consider the gen- 
eral case that a depends on time - since all the quan- 
tities in Eq. (|l^ ) depend on time -, but find that a in- 
creases just negligibly in our models, thus recovering the 
KB90 assumption that the velocity of the shell is almost 
constant during the momentum-driven phase. An exam- 
ple is shown in Fig. ^ where Vagb — 15.0 km s^^ and 
Vs ~ 16.3 km s~^, in very good agreement with the re- 
sults of hydrodynamical calculations by Kifonidis (1996). 
The plot also shows the predicted speed-up of the fast 
wind (with velocity v according to the prescriptions in 
Sect. 2.3.1) and the escape velocity Wcsc during the mo- 



mentum driven phase. 

The momentum- driven phase ends when the hot 
shocked gas is able to fill the interior in a time shorter 
than the cooling time. Following KB90, the change-over 
time from the momentum-driven phase to the energy- 
driven phase, crucially depends on the velocity v of the 
fast wind: 

A/agb qa 



ij^ — 5.23 



(1 + V^) 1 



(15) 



quantities are in c.g.s. units. We obtain that Eq. (|15|) is 
fulfilled typically for v ^ 100 — 150 km s^^, in agreement 
with the findings of KB90. Clearly, is quite sensitive 
to the evolutionary speed of the central star during the 
initial post-AGB stages, as illustrated in Fig. ^ for two 
choices of the transition time. 

Finally, it is worth remarking that in all our models, 
during the momentum-driven phase, the shell remains es- 
sentially neutral, the ionisation of H possibly starting to- 
wards the end of this phase and involving an almost neg- 
ligible part of the nebula. 

2.4.2. The energy-driven phase 

The dynamical evolution of the shell at t > is mainly 
controlled by the kinetic energy input from the fast wind, 
which can be substantially thermalised as radiative cool- 
ing is inefficient. A hot bubble develops, and the physi- 
cal conditions of the interior gas are those produced by a 
strong adiabatic shock. The front of the inner shock sits 
very close to the central star (i.e. i?i <C i?s)- During this 
phase the thermal pressure inside the hot bubble is the 
driving force acting on the shell. 

In our model, as soon as t = we switch to the dynam- 
ical description developed by VK85. From the conditions 
of conservation of mass, momentum, and energy applied 
to the system consisting of the hot bubble and the shell, 
one finally ends up with a set of coupled, non-linear, first- 
order differential equations: 

At 



dt 

dt 

dP 

dT 



= A^AGB 



Vs 



- 1 



A^RlP - A/agb 



(16) 
(17) 

(18) 

MRl-Rl) ^^^^ 

Radiative cooling in the shocked region is neglected, which 
is a reasonable approximation during the energy-driven 
phase. The numerical integration of Eqs. (|l^) - (|l9| ) is 



V Vagb 

1 

Ms 

(l/2)mw2 - 107rPsVsi?i 



(Vs - Vagb? 
Vagb 
67rPsVii?i2 
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performed explicitly, with a integration step much shorter 
(i.e. 1/100) than the actual evolutionary time step. This 
latter is chosen so that the maximum variation cither in 
the effective temperature or in the luminosity of the cen- 
tral stars cannot be larger than 0.01 dex. The initial con- 
ditions for integration are determined by the values of the 
variables at the end of the momentum-driven phase. 

It is worth remarking that in our model, both Magb 
and Vagb are not constant as assumed in many past 
works, but are given as a function of time according to the 



synthetic calculations of the AGB phase (see Sect. 2.2). 
This means that, at each instant, A/agb and Vagb repre- 
sent the mass-loss rate and the expansion velocity of the 
wind layer previously ejected during the AGB, which is 
now being swept by the shocked shell. 

In Eq. ( p^ ) the radius (i?i) and velocity (Vi) of the 
inner shock also enter. It was shown by hydrodynamical 
calculations (e.g. Marten & Schonberner 1991) that i?i 
rapidly decreases from its initial value (^ i?s, as in the 
momentum-driven phase) to a value small compared to 
Rs- Therefore, it is reasonable to adopt for Ri at each 
time step the asymptotic solution as derived by VK85: 



3l 

Rs 



Srhv 



4MagbVagb VVagb 



At each instant, the velocity Vi 
Eq. m. 



(20) 



di?i/dt, is derived from 



2.4.3. Dynamical interaction of AGB shells 

As pointed out by hydrodynamical simulations (e.g. Frank 
et al. 1990, 1994; Steffen & Schonberner 2000; Corradi et 
al. 2000), beyond the primary shocked shell blown by the 
fast wind, dynamical interactions also occur in the "undis- 
turbed" AGB wind, which may account for a number of 
complex morphological and emission features, i.e. crowns, 
edges, multiple shells, extended halos, etc. It is obvious 
that a similar detailed analysis is out of reach for our ana- 
lytical model, nevertheless we can at least attempt to give 
a zero-order description of AGB wind interactions. 

According to the adopted formalism for the AGB wind 
(VW93), both M and Vagb increase with time until the 
natural onset of the superwind regime. In this scenario, it 
turns out that interactions between TP-AGB ejecta with 
different expansion velocities may be possible, even if they 
are not yet reached by the outer shock front produced by 
the fast wind. 

Let us consider three consecutive AGB shells, with 
masses Mi, M2, M3, that are expanding with constant 
radial velocities at a given instant t. For the sake of sim- 
plicity, we denote with {Vi\ i = 1,4) the velocities at the 
shell boundaries of radii {Ri; i — 1,4), so that the i*'^- 
shell is moving with velocity Vi. The situation is sketched 
in Fig. H Suppose also that V2 > V3. As a consequence, it 
may happen that at a certain time t+dt during the expan- 
sion, shell 2 catches up shell 3. This would cause a sort of 
overlap of gas shells with different dynamical properties. 









1 


2 


3 



3 2 

m; m. 



radial coordinate 

Fig. 5. Schematic representation of the possible interac- 
tion between expanding shells ejected during the AGB. 
See text for a complete description 

and hence a change in the mass and velocity distributions 
along the radial coordinate. 

To describe the final result (at time t+dt) of this inter- 
action we adopt the following simple scheme. We denote 
with 



(21) 



the new radial location of the shell boundaries (at increas- 
ing distance; note the exchange between indices 2 and 3), 
and assume that V{ = V\ and = V4 (i.e. at these radii 
the expansion motion is undisturbed). Then, the problem 
is to determine the five unknowns {M[, Mj, M3, V2, and 
Kj') which characterise at time t + dt the shells delimited 




by radii i?^-, as given by Eqs. (|2l|). From imposing the mass 
conservation of the system composed by the three original 
shells, we get the new masses: 




Ml 

M2 
M3 



AMi 
AMi 
AM3 



/ 3 



AM, = 



r?/ 3 7-)/ > 



R'a 



R'. 



Ml 



M3 



(22) 



(23) 



whereas from the conditions of momentum and energy 
conservation: 

MiVi + M2V2 + M3V3 = M[Vl + M2V2 + M3T/3 (24) 
AhVi^ + M2V2^ + MsVs^ = M[vf + M^vf + Ml^vf 



p. Marigo et al.: Evolution of Planetary Nebulae I. An improved synthetic model 



9 



-15 



Q. -20 

QO 
O 

-25 



-15 
C5, -20 

O 

-25 

-15 

Q. -20 
o 

-25 




t = 60 



t= 15041 




t=30142 




10 L 

- Il-burners 



2 4 6 

radius (pc) 

Fig. 6. Evolution of the radial density (in gr cm"'^) struc- 
ture of the circumsteUar ejecta (model 4; see Table H) at 
selected ages (indicated in yr) during the post-AGB phase. 
t = would correspond to the moment in which the cen- 
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we obtain the new velocities V2 and (for which we omit 
to present the formulas due to their trivial derivation). 

This produces evident perturbations in the density 
pattern that would be otherwise expected for a station- 
ary gas flow. A possible consequence of the shell overlap 
sketched in Fig. || is the appearance of a density enhance- 
ment within the new shell 2 (with mass M2+ AMi+ AM2) ■ 

According to the above prescriptions we find that this 
kind of dynamical events begins to affect the circumsteUar 
ejecta in the last stages of the AGB evolution (i.e. after 
the onset of the super- wind regime) , and usually involves 
shells ejected during the transition from the slow AGB 
wind to the superwind. This feature is in agreement with 
the hydrodynamical findings by Steffen & Schonberner 
(2000). 

An example is shown in Fig. ||. We can notice that 
non-negligible density enhancements show up in the radial 
density structure of the model. These density spikes are 
rather close to the primary shocked shell, while the long 
tail of slow low-density shells remains unperturbed. 

We are completely aware that such approach is an ex- 
tremely rough description of the real situation, which can 
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[ ■ Clegg & Middlemass 1987 



He— burners 




log T.u 

Fig. 7. Available data grid for the ionising rates of H and 
He^ in the log Toff — log g plane according to model atmo- 
spheres of Clegg & Middlemass (1987; filled squares) and 
Kurucz (1994; empty squares). The sets of H- and He- 
burning post-AGB tracks from Vassiliadis & Wood (1994) 
are superimposed 



be handled correctly only with the aid of hydrodynamics. 
Nevertheless, it may be interesting to get a first hint of 
likely interaction events (i.e. when and where they might 
occur) in the unperturbed AGB ejecta surrounding the 
primary shocked nebula. Furthermore, we remark that the 
dynamical scheme presented in this section does not affect 
any of the relevant predictions, as long as the PN is as- 
similated to the primary shocked s hell, as in the analytical 
interacting- winds model (see Sect. 2.4). 



2.5. Photoionisation rates for H and He^ 

As the central star evolves at increasing effective temper- 
atures along the horizontal part of the H-R track, it starts 
emitting UV photons that can ionise the surrounding neb- 
ula. 

In this study we account for the photoionisation of H 
and IIe+, which occurs when a significant fraction of the 
emergent photons has a wavelength shorter than 911 A 
(for Teff ^ 30 000 K) and 228 A (for T^g ^ 60 000 K), re- 
spectively, lonisation fluxes (in cm^^ s~^) are derived from 
atmosphere models, namely: the LTE models calculated 
by Kurucz (1994) in the ranges: 10 000 < Toff < 40 000 K, 
2.0 < logg < 5.0; and the non-LTE models by Clegg 
& Middlemass (1987) in the ranges: 40 000 K < T^s < 
180 000 K, 4.0 < logg < 8.0. In the adopted model atmo- 
spheres the abundance ratio IIe/H= 0.1 (by number) is 
assumed. 
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Fig. 8. lonisation rates as a function of time for a 0.6 Mq 
post-AGB central star (model 4 in Table ||) . See text for 
more details 

In Fig. 1^ the adopted grids of model atmospheres are 
shown in the log Tcff — log g plane, together with the sets of 
post-AGB tracks for both H- and He-burner central stars. 
For each grid point, ionising fluxes are obtained integrat- 
ing the synthetic spectra 

Q = ^A{hv) (25) 

with hvi — 13.6 eV and hv2 — 54.4 eV for the H con- 
tinuum, and hvi — 54.4 eV and hv2 = oo for the He"'' 
continuum, thus neglecting absorption by neutral helium 
in the nebula. These calculations yield Q{BP) and Q(Hc^) 
as a function of both effective temperature (in the range 
Teff : 10 — 180 X 10'^ K) and surface gravity (in the range 
\ogg : 2.0 - 8.0) of the hot star. 

Within the g — Tcff range spanned by the model at- 
mospheres, a bi-linear interpolation in log Tcff and log 5 
over this data grid is then adopted to evaluate the current 
ionisation fluxes along the post-AGB tracks. For possible 
points not covered by the grid, we extrapolate the model 
data in \ogg for Teff < 180 000 K, and we assume a black- 
body spectral energy distribution for Toff > 180 000 K 
(temperatures attained only by H-burning central stars 
with M Z 0.7Mq; see Fig. @). 

An example is shown in Fig. ^, where the ionisation 
rates for H and He+ (i.e. ^ttRIQ{YL) and 47ri?2Q(He+) in 



s^^, where i?* is the stellar radius) are calculated as a 
function of time for an evolving central star with a mass 
of about 0.6 Mq. Predictions from Kurucz and Clegg & 
Middlemass model atmospheres are compared to black- 
body integrations. 

The temporal behaviour of the ionisation rates is con- 
trolled by the changes in Tcff and L, and the evolution- 
ary speed of the central stars. In Fig. ^ we can see that 
the photoionisation rates significantly increase as the star 
heats up at nearly constant luminosity, reach a peak (in 
proximity of the knee of the H-R track, at the maximum 
of Tcff in the top panel) and then quickly fall down when 
the luminosity drops. Subsequently, the slower decrease 
of the ionisation rates reflects the drastic deceleration of 
the evolutionary speed. We can also notice that in general 
(5(11) > (3(IIe+), and the peak duration is larger for (5(H) 
compared to Q{]le^). 

As far as the ionisation of H is concerned, we notice 
that the blackbody assumption notably overestimates the 
continuum flux for Tcff < 40 000 K, whereas it approxi- 
mates well the atmosphere models at higher temperatures. 

Much larger differences exist between the He+ ionisa- 
tion fluxes predicted according to the blackbody and stel- 
lar atmospheres models. These discrepancies are due to 
important non-LTE effects that significantly alter the rel- 
ative number of neutral, single, and doubly ionised He ex- 
pected under LTE conditions (Clegg & Middlemass 1987). 

As shown in the bottom panel of Fig. ||, for most of the 
post-AGB evolution the adoption of a blackbody spectrum 
leads also to overestimate the Hc+ continuum, except for 
the emission peak corresponding to the highest effective 
temperatures, before the decline of the stellar luminosity 
(see top panel). Such remarkable excess of UV photons 
may be one of the main factors at the origin of the so- 
called Zanstra discrepancy (refer to Gruenwald & Viegas 
2000 for a recent analysis; see also Sect. 4.4). 



2.6. Ionised mass 

The ionisation structure of the nebula is solved with the 
aid of simple but reasonable prescriptions. Once the ionis- 
ing rates emitted by the central star are known, the ionised 
mass can be singled out, provided that the current density 
stratification of the circumstellar gas (i.e. shocked shell 
and unperturbed AGB shells) is determined. The density 
of the ionised part of the shell is estimated from the equa- 
tion of state of a perfect gas, assuming a constant electron 
temperature To (typically 10 000 K). 

Then, the ionised masses for both H and He+ are de- 
rived assuming equilibrium between ionisation and recom- 
bination events. It means, that within a certain volume, 
limited by the Stromgren radius i?ion, the rate of recom- 
binations to all levels but the ground state, is equal to the 
number of ionising photons from the central star. This 
condition translates into 



/■i?,on(H ' ; 

Q(H°) = 47r / aB{ii°yN^N{Il+)dr 
Jr.* 



(26) 
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for neutral hydrogen and 
Q(He+) = An 



aB(Hc+)r2iVoiV(He++)dr (27) 



for single-ionised helium. It is clear that these equations 
may be fulfilled as long as i?ion < -Knob, where i?ncb is the 
outermost radius of the gaseous nebula. The quantities ae 
are the effective recombination coefficients, expressing the 
total recombination rates (cm^ s~^) to all excited levels. 
They depend weakly on the electron temperature To and 
can be expressed as 



aB(H°) 
aB(He+) 



1.5210-12 (r,/io4)-o-7i^ 



(28) 



according to the fitting formulas presented by Oliva & 
Panagia (1983). 

We assume that the ionisation radius for single- 
ionised helium coincides with that of ionised hydrogen, 
i.e. i?ion(He^) = i?ion(H+), which is a good approxima- 
tion for T ^ 4 X 10^ K (Osterbrock 1974). Then, denoting 
by A^(H) and A^(He) the total number densities (neutral 
atoms -t- ions; in cm"-^) of hydrogen and hehum, and by 
A'e the number density of free electrons, we adopt the fol- 
lowing scheme for the ionisation structure (with obvious 
meaning of the particle number densities TV): 



iV(H") = 7V(He°) = iV(He+) = 
iV(H+) = N{B) 
iV(He++) = Ar(He) 
iVo = iV(H) + 2 7V(He) 

for r < i?ion(He++); 

iV(HO) = TV (He") = iV(He++) = 
iV(H+) = 7V(H) 
iV(He+) = iV(He) 
iVc = N{K) + AT (He) 

for i?ion(Hc++) < r < i?io„(H+), and 

' iV(H°) = N{K) 
iV(H+) = Ar(He+) = iV(Hc++) = 
7V(He°) = 7V(He) 
iVe-0 



(29) 



(30) 



(31) 



for r > i?ion(H+). 

2.7. Shell thickening due to ionisation 

An implicit assumption of the interacting-winds model 
(outlined in Sects. 2.4.1 and 2.4.2) is that the shocked shell 
is thin compared to its radial extent, i.e. ARg/Rs ^ 1- 
However, observations indicate that ring-like nebulae ex- 
hibit a large range of relative thickness (from about 0.2 
to 0.8, with the peak of the distribution at ~ 0.4 — 0.5), 
and even homogeneous nebulae lacking a central hole ex- 
ist (Zhang & Kwok 1998). Moreover, the expansion ve- 
locities of many young PNe can be notably high, up to 
~ 50kms-i (e.g. Weinberger 1989), which is something 



not predicted in the framework of the simple interacting- 
winds model. 

Hydrodynamical models of planetary nebulae (Marten 
& Schonberner 1991) have pointed out that these observa- 
tions can be explained considering the dynamical effects 
caused by the ionisation of the nebula. In fact, thanks to 
the growing number of free electrons the thermal pressure 
increases inside the ionised nebula, which then tends to 
expand both inward and outward. As a consequence, the 
inner and outer nebular rims are decelerated and acceler- 
ated, respectively. 

Such dynamical effect is included in this model follow- 
ing the analytical treatment suggested by Kahn (1983). 
The basic steps are shortly recalled below. In practice, 
the effect of ionisation is treated as a linear perturbation 
of the equations of momentum and energy conservation 
(Eqs. 18-19|), by introducing a term APion for the ex- 



tra pressure. Once the unperturbed solution of the set of 
Eqs. ([l^ ) - ( [l9| ) is singled out (in the thin-shell approxi- 
mation) , proper first-order corrections are then applied to 
the radial coordinate Rs and velocity Vs of the contact 
discontinuity according to the following scheme. 

Let us assume that, as a consequence of the pressure 
contribution due to ionisation, the neutral shell expands 
outward to a radius Rs + Ai?a (with Ai?a > 0), and the 
contact discontinuity expands inward down to a radius 
Rs + Ai?b (with Ai?b < 0), where Rs refers to the unper- 
turbed value. Then, the unperturbed energy equation 



^(2.i^iP) 

becomes 
d 



dt 
1 



—mv 
2 



AttRIP 



dRs 
dt 



[2^(i?s + Ai?b)^(P + APio„)] = 
- 47r(Ps + ARkfiP + APio 



d(Ps + APb) 
dt 



(32) 



(33) 



An analogous treatment is applied to the unperturbed mo- 
mentum equation. The relationship between Pa and Pb is 
given by equation 31 in Kahn (1983) to whom we refer for 
all details. Expliciting the time derivatives and neglecting 
the second-order terms, one eventually derives the radius 
corrections (APa, APb), the velocity corrections (APa, 
APb), and the pressure corrections APjon, APjon as a 
function of the unperturbed quantities, winds parameters 
and ionised mass. 

For the sake of clarity, we can now summarise our pre- 
scriptions. As the ionisation front breaks into the shocked 
shell, this latter will be characterised by three relevant 
radii, namely: 



Rg — Rs + APb : actual contact discontinuity 



(34) 



Pion = Rs + APa : ionisation front 

1/3 



Pont — 



Pir 



3Mn 



(35) 

: outer radius of the shell(36) 



The outer radius Pout of the neutral part of the shell - 
with mass M^out = Ms — Afion - is calculated assuming a 
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Fig. 9. Predicted properties of the 
shocked nebula as a function of 
time, accounting for the dynam- 
ical effect resulting from ionisa- 
tion (model 4 in Table |). Top- 
left panel: Temporal behaviour 
of the mass of the shocked shell 
Afg, and ionised masses for H"*" 
and He+~''. The total mass Mmax 
of the circumstellar AGB ejecta 
is also shown (dotted line). Top- 
right panel: Comparison between 
the rate of injection of mechanical 
energy from the fast wind (m?;^/2), 
and the rate of the work done by 
the pressure in expanding the shell 
(— i^w, — -E-w^ see text). Bottom- 
left panel: Relative thickness 
of the ionised part of the shocked 
shell. The prediction under the 
thin-shell approximation is com- 
pared to the case in which the 
shell thickening due to ionisation 
is included. Bottom-right panel: 
Velocities at the shell boundaries: 
of the contact discontinuity un- 
der the thin-shell approximation 
(dotted line); Vg of the actual con- 
tact discontinuity, and T^on of the 
ionisation front (solid lines) when 
the shell-thickening corrections are 
included. See text for more details 



uniform average density Pncut, which is derived from the 
equation of state of perfect gas, assuming a temperature 
of 100 K. In this case the isothermal speed of sound is 
typically of 1 km s~^. 

In general, we have i?g < i?ion < -Rout, the first and 
second equalities corresponding to the cases of complete 
shell neutrality (Mneut = Afg) and ionisation (Mnout = 0), 
respectively. 

Moreover, it follows that the expansion of the ionised 
part of the shell produces a deceleration of the actual con- 
tact discontinuity that will expand with a velocity 

Vg* = + Ai?b (with A_Rb < 0) , (37) 

and an acceleration of the ionisation front that will move 
outward with a velocity 

Mon = + Ai?a (with Ai?a > 0) (38) 

It should to be noticed that such a first-order per- 
turbation approach is strictly correct only as long as 
Ai?b .Rs, i-e. during the early stages of the expansion 
when a small fraction of the shell is ionised. Moreover, with 
the above prescriptions we implicitly suppose that there 
is a pressure equilibrium and homogeneity in the ionised 



part of the shell. This assumption is acceptable consider- 
ing that a rarefaction wave, generated by the initial inward 
motion of the contact discontinuity, is expected to rapidly 
homogenise the gas (Breitschwerdt & Kahn 1990). It is 
also clear that the adopted treatment is a very simpli- 
fied approximation of the real phenomenon. Nevertheless, 
comparing our results with both the predictions of hydro- 
dynamical models and observed data, we generally find a 
good agreement (see Figs. || and |T7| ). 

An example of the application of the shell-thickening 
approximation is shown in Fig. ^, that refers to a model 
nebula surrounding a central star with mass of ^ 0.6 Mq, 
and initial solar metallicity. 

The relative thickness of the ionised shell 

bR = ^'""'^^ (39) 

is shown as a function of time (solid line, bottom-left 
panel), and compared to that expected according to the 
thin-shell approximation (dashed line). The latter curve 
is derived by integrating Eqs. (|l^) - (|l^), and assuming 
that the gas density in the ionised shocked shell is a fac- 
tor k = h times larger than the density of the unperturbed 
AGB wind, i.e. pion = ^AfAGB/(47ri?gVs), as suggested by 
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VK85.[] In this case the thickness slowly increases with 
time, being mostly below 0.1. 

On the contrary, when the pressure contribution due 
to ionisation is taken into account, the shocked shell pro- 
gressively thickens as the ionisation front breaks into it, 
reaching a maximum of SR ^ 0.7 when the whole shell is 
ionised. Actually, at this stage, the ionisation front out- 
passes the outer rim of the shocked shell and embraces 
the undisturbed low-density AGB wind completely (see 
top- left panel of Fig. |9[). 

After the maximum of SR, its subsequent decline and 
further increase can be understood considering the compe- 
tition between the rates of i) the mechanical energy input 
from the fast wind - which tends to compress the shell 
from the back (term - and ii) the work done by 

the pressure to expand the shell. The latter is expressed 
by the term —Ew = 47r/'i?gVs for the unperturbed case, 
and the equivalent term —E^ — Att{P + APion)(i?s + 
Ai?b)^(Vs -|- Ai?b) when the first-order corrections are ap- 
plied. As shown in Fig. || (top- right panel), as long as the 
first term dominates, the shell thickness decreases; at later 
times it starts to increase again, as the kinetic energy in- 
put of the fast-wind drops and the adiabatic expansion 
prevails. 

Similar considerations can be applied to explain the 
temporal behaviour of the velocity at the inner and outer 
rims of the ionised nebula (bottom-right panel). The in- 
crease of the pressure inside of the ionised shell determines 
the initial strong acceleration of the ionisation front and 
the deceleration of the contact discontinuity (which may 
even attain negative values corresponding to a short-lived 
inward motion). The subsequent acceleration of both rims 
(the shell is now fully ionised) is driven by the increasing 
pressure of the hot bubble. Finally, when the strength of 
the fast wind weakens, the contact discontinuity and the 
outer rim of the nebula are decelerated and accelerated, 
respectively, while the ionised gas tends to expand both 
inward and outward. We can note that the velocities at- 
tained at the same time by the shell rims can be quite 
different, and cover the wide range of observed velocities 
in PNe, typically from few km s~^ to 30 - 40 km s~^. 

The reliability of these predictions can be checked by 
comparison with the results of hydrodynamical calcula- 
tions from Marten & Schonberner (1991). As far as the 
shell thickening is concerned, it is encouraging to note 
that the temporal behaviour of SR shown in Fig. ^ closely 
resembles the predictions of their numerical simulations 
(see figure 10 in Marten & Schonberner). The evolution of 
the rim velocities - i.e. the acceleration of the ionisation 
front and deceleration of the contact discontinuity during 
the early stages of the ionisation phase - is also recovered 
by our analytical model. 



2.8. Optical properties of the nebulae 



With the aid of the prescriptions given in Sect. 2.6, we 
are able to determine not only the ionised masses for H+, 
He+, and He++ at each time step during the post- AGB 
evolution, but wc can also naturally discriminate between 
optically thick and optically thin configurations, depend- 
ing on whether the continuum photons (relative to a given 
ion) are partially or totally absorbed by radiative recombi- 
nation transitions in the nebula. To this aim, we calculate 
at each time during the evolution the absorbing factor (as 
defined by Mendez et al. 1993) 



N, 



abs 



(40) 



^ fc = 5 is indicated by VK85 as a reasonable compromise 
between the high density contrast generated by a isothermal 
shock and the density decrease in the ionised nebula 



giving the fraction of the continuum emitted photons that 
is actually absorbed by the ion species under considera- 
tion. In practice, it can be expressed as the ratio between 
the ionisation rate Q and the recombination rate within 
the Stromgren radius, i?ion (see Eqs. |2^). It follows 
that /i = 1 as long as i?ion < -Rncb, which corresponds to 
the optically thick case. On the other hand, if no radial co- 
ordinate, up to i?nob, can be found where the equilibrium 
condition is satisfied - so that i?ion cannot be formally de- 
fined -, we have /i < 1, which corresponds to the optically 
thin case. Figure |l^ illustrates how the nebular optical 
properties in both H and He+ continua are predicted to 
vary as a function of time as the central star evolve on its 
post-AGB track. 

As far as the 11+ continuum is concerned, we can see 
that when the central star becomes hot enough to ionise 
the H atoms, the nebula is expected to be initially opti- 
cally thick. The transition to the optically-thin configura- 
tion occurs as the effective temperature keeps increasing 
along the horizontal part of the track and the number 
of the emitted ionising photons is growing. At this stage 
the whole nebula is ionised. In proximity to the knee of 
the track (the point of maximum temperature), fj, attains 
its minimum value, implying that a large number of pho- 
tons are escaping from the ionised nebula. Subsequently, 
the luminosity starts to drop, the photoionisation rate de- 
creases (hence ^ increases), and the nebula becomes opti- 
cally thick again. 

Later on, a second conversion to the optically-thin con- 
figuration occurs, being essentially the result of the dy- 
namical expansion of the nebula. In fact, when the gas 
density gets low enough, recombination in the whole neb- 
ula becomes less efficient (in terms of rate of events) than 
photoionisation. 

The optical properties in the He++ show notable dif- 
ferences compared to the case of H"'". First, helium starts 
to be doubly ionised at quite high effective temperatures 
(typically Toff ~ 60 000 K). Second, the optically-thin con- 
figuration is very short-lived, and restricted to the knee of 
the track. This reflects the fact that the photoionisation 
rate for Hc+ is systematically lower than that for H (see 
Fig. H), and presents a much steeper decrease with effec- 
tive temperature and luminosity. 
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Fig. 10. Expected optically thick 
and thin configurations during the 
evolution of a PN ionised by a 
central star with a mass of 
0.6Mq and initial solar composition 
(model 4 in Table ||). Left panels: 
absorbing factor ^ (given by Eq. |o|) 
as a function of time for both H+ 
and He"*""*" ions (top and bottom 
panels, respectively). Time is set to 
zero at the end of the AGB evolu- 
tion. Right panels: Predicted op- 
tical appearance of the nebula as 
the central star evolves along the 
post-AGB track on the H-R dia- 
gram. The adopted notation is as 
follows: dotted line as long as the 
photoionisation rate Q < lO'*" s^^; 
solid lines for optical thick (^ = 1) 
or thin (^ < 1) configurations as 
indicated. The horizontal mark on 
the tracks corresponds to the evo- 
lutionary age of 4 X 10* yr, which 
corresponds to the limit of the in- 
terval shown in the left plots 



Finally, it is worth remarking that in our model ionisa- 
tion and optical properties of PNe are predicted according 
to a coherent coupling between the evolution of the central 
star and the dynamics of the nebula. This represents an 
important improvement upon past synthetic approaches, 
that in many cases lack any dynamical description of the 
nebulae (e.g. Stanghellini & Renzini 2000) and/or adopt 
fj, as free parameter (e.g. Mendez & Soffner 1997). 

2.9. Intensities of few relevant optical emission lines 

We aim at predicting the nebular luminosities at the 
wavelengths of a few important lines typically observed 
in PNe, namely: the recombination lines H/3 A4861 and 
He II A4686, and the nebular forbidden line [O ill] A5007. 

As a first approach, the intensities of these lines are 
calculated by standard procedures (as described below), 
under simple but reasonable assumptions for the physical 
conditions of the emitting regions. Specifically, we assume 
that the temperature Tc and density A'e of free electrons 
are constant. In particular, is predicted according to 
the ionisation structure described by Eqs. ( |2^ ) - (^l|), 
whereas To is considered a free input parameter, which 
is chosen in the typical range indicated by observations, 
that is 8 000 K ^ Te ^ 15 000 K. 

These assumptions may be considered completely sat- 
isfactory as far as recombinations lines are concerned. 



as they depend mostly on atomic parameters and only 
weakly on To. The case of forbidden lines is more delicate, 
due to their high sensitiveness the electron temperature 
which cannot be predicted but has to be assumed by the 
model. It is clear that an improved description of the line 
emission from PNe can be attained only with the aid of 
a photoionisation code. However, our choice should repre- 
sent a reasonable first step, and should be adequate to the 
purpose of exploring the sensitiveness of PNe properties 
to various possible physical conditions. 

2.9.1. The recombination lines: H/3 and HeiiA4686 

The evaluation of the total luminosity Lx emitted by the 
ionised nebula in a particular recombination transition (of 
wavelength A) requires the calculation of 

Ta = 47r / Jxir^dr (41) 
Jr* 

where the integral is performed over the Stromgren volume 
of the ion under consideration, with 

Jx = a{nn') Ni^n hv (42) 

being the volume emission (erg cm~'^ s"^), and ainn') 
denoting the recombination coefficient of the transition 
from state n to state n' . 
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This latter parameter is calculated adopting the fit- 
ting relations presented by Pequignot et al. (1991), as a 
function of Tc and ionic charge, for both the recombina- 
tion H/3 line of H at A = 4861 A [aH+(4,2)], and the 
recombination line of He ii at A = 4686 A [aHc++ (4, 3)]. 

All other quantities in the above equations are known 
once the ionisation structure of the nebula is determined 



following the prescriptions described in Sect. 2.6 



2.9.2. The forbidden line: [0 iii] A5007 

The intensity of this nebular line - corresponding to the 
transition from the coUisionally-excited metastable level 
j — to the level i — - is calculated from 



r(o^ 



r(0) 



no), 

Y(H) " 



l + y(He)/y(H) 



(43) 



(expressed in erg cm"'^ s"^), where Aji = 0.021 s^^ is 
the radiative transition probability; Nj(Te, N^) is the rela- 
tive population of level j; Y{0^~^)/Y{0) denotes the ionic 
abundance of 0++ relative to the total oxygen abundance 
(in all excitation and ionisation states); Y{0)/Y{}i) refers 
to abundance of oxygen (in all states) with respect to 
hydrogen; Y{R) and ^(He) are the chemical abundances 
(by number, in mole gr~^) of hydrogen and helium; and 
n = N(H) -|-7V(He) corresponds to the number density (in 
cm~^) of H and He. 

The most critical quantity here is Nj{Te, N^), as it de- 
pends on both the electron temperature and density. It is 
determined using the detailed calculations of atomic tran- 
sitions performed by Kafatos & Lynch (1980), interpolat- 
ing in Te and Nc the data of their Table 3E. All other 
quantities (i.e. elemental abundances and densities) are 
predicted by our model. 

Equation (|4^) must be integrated over volume to ob- 
tain the emitted luminosity in the [O ill] A5007 line. We 
assume that both 0++ and He"*" ions share the same 
Stromgren volumes, based on the fact the ionisation po- 
tentials of 0++ (54.9 eV) and He+ (54.4 eV) are almost 
identical. In other words, we assume that the emission of 
the [O III] A5007 line can only originate from regions where 
He is single-ionised, disappearing as soon as He becomes 
double-ionised. Moreover, we consider that within that re- 
gion Y{0^~^) = Y{0), which is also a good approximation 
(Osterbrock 1974). 

2.10. Summary of model parameters 

Before analysing the expected dynamical and ionisation 
properties of PNe, it is worth summarising the parame- 
ters/prescriptions of the model that are can be varied in 
our calculations, namely: 

— initial metallicity of the stellar progenitor, Zi; 

— stellar mass at the onset of the TP-AGB phase, A^itp; 

— mass of the central star of the PN model, Mcspn; 

— H- or He-burning post-AGB track; 

— the transition time, ttr] 



— expansion velocity of the stellar wind during the last 
stages of the AGB evolution, V^g; 

— electron temperature of the ionised shell. To. 

It is worth remarking that our default set of prescrip- 
tions include i) the ionisation fluxes from model atmo- 
spheres, and ii) the thick-shell correction scheme. For pur- 
pose of comparison, we may easily relax these assumptions 
adopting the blackbody energy distribution and the thin- 
shell approximation. This is the case for few examples dis- 
cussed in the text. 

Calculated PNe models are listed in Table |[ so that 
each reference number Nmod corresponds to a particular 
set of parameters. It is clear that some of these quantities 
are intrinsic properties of the models (e.g. Zi, Mcspn), 
whereas others are genuine free parameters (e.g. ttr, I^agb; 
Tc). The aim is to define suitable ranges for these parame- 
ters, on the basis of observations of Galactic PNe. For the 
other model ingredients, not explicitly mentioned here, we 
refer to the prescriptions already discussed in Sect. |[ 

3. Overview of model predictions 

3.1. Summary plots 

The several panels in Fig. |ll] illustrate the time evolution - 
during a period of 30 000 yr - of a PN model corresponding 
to a H-burning central star with a mass of about 0.6 Mq. 
The selected model (number 4 in Table ||) may be consid- 
ered typical of Galactic PNe. The run of several physical 
quantities is shown simultaneously, then offering a useful 
guide to our further discussion. In short, going from top 
to bottom panels, we present: 

1. the evolution of luminosity and temperature for the 



central star (left; Sect. 2.3), and the corresponding ion- 
isation rates (right; Sect. 2.5); 

the masses of the shocked shell and ionised re- 
gions (left), and the corresponding fi factors (right; 
Sect. |?|); 

the radii (left) and velocities (right) of the shell bound- 
aries ( incl uding the radial location of the inner shock) 
(Sect. 13); 

the line ratios (left) and total luminosities (right) in 
the main emission lines (Sect. UM. 



Similar plots are available for all our models, and 
allow a better interpretation of the different evolution- 
ary features that may be present. A very interesting 
comparison is provided, for instance, between the H- 
burning track of Fig. ^l], and a model computed with 
the same set of parameters, but using a He-burning 
track (number 9 in Table |^). The latter is presented 
in Fig. for a total period of 55 000 yr. As can be 
readily noticed, the He-burning track is characterised by 
a smaller luminosity and slower evolution than the H- 
burning one, and this determines a substantially differ- 
ent evolution of the nebular properties. For instance, the 
nebula shown in Fig. [13 is nearly always optically thin 
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Table 2. Specification of the basic parameters adopted to calculate PN models 



AT 

"mod 


AT 


JVJCSPN 


xi-/ xie- 


T /max 
VAGB 




T' 


notes 








burner 


t^Kin s ) 


[W yrj 






1 


1.7 


0.5989 


H 


15 


1.5 


10.0 


c 


2 


1.7 


0.5989 


H 


f{Z) 


1.5 


10.0 


c 


3 


1.7 


0.5989 


H 


fiL)* 


1.5 


10.0 


c 


4 


1.7 


0.5989 


H 


.f{L) 


1.5 


10.0 


A, B, C, D, E 


5 


1.7 


0.5989 


H 


15 


0.5 


10.0 


C, D 


6 


1.7 


0.5989 


H 


f{L) 


1.5 


10.0 


thin-shell approx. 


7 


1.7 


0.5989 


H 


f{L) 


1.5 


10.0 


blackbody 


8 


1.7 


0.5989 


H 


.f{L) 


8.0 


10.0 


D 


9 


1.7 


0.5989 


He 


f{L) 


1.5 


10.0 


B, D 


10 


1.7 


0.5989 


He 


f{L) 


5.0 


10.0 


D 


11 


1.7 


0.5989 


H 


f{L) 


1.5 


8.0 


E 


12 


1.7 


0.5989 


H 


f{L) 


1.5 


12.0 


E 


13 


1.7 


0.5989 


H 


fiL) 


1.5 


13.0 


E 


14 


1.0 


0.5577 


H 


f{L) 


10.0 


10.0 


C 


15 


1.0 


0.5577 


H 


f{Z) 


10.0 


10.0 


C 


16 


1.2 


0.5706 


H 


f{L) 


1.5 


10.0 


A, D 


17 


1.2 


0.5706 


H 


f{L) 


8.0 


10.0 


D 


18 


1.5 


0.5868 


H 


fiL) 


1.5 


10.0 


A 


19 


1.5 


0.5868 


He 


f{L) 


10.0 


10.0 




20 


1.8 


0.6120 


H 


.f{L) 


1.5 


10.0 


A 

A 


21 


2.0 


0.6296 


H 


f{L) 


1.5 


10.0 


A, B ,D 


22 


2.0 


0.6296 


H 


f{L) 


0.5 


10.0 


D 


23 


2.0 


0.6296 


He 


f{L) 


1.5 


10.0 


B 


24 


2.2 


0.6529 


H 


fiL) 


1.5 


10.0 


A, D 


25 


2.2 


0.6529 


H 


fiL) 


0.5 


10.0 


D 


26 


2.5 


0.6865 


H 


fiL) 


1.5 


10.0 


A, D 


27 


2.5 


0.6865 


H 


f{L) 


0.5 


10.0 


D 


28 


2.5 


0.6865 


H 


f{L) 


0.3 


10.0 


D 


29 


3.0 


0.7867 


H 


f{L) 


0.5 


10.0 





Sequences of models with the same set of model parameters, but varying 

A: McsPN (4-16-18-20-21-24-26), (9-16); 

B: H-/He-burning post-AGB tracks (4-9), (8-16); 

C: V^gI (1-2-3-4-5), (14,15); 

D: Ttr (4-8), (9-10), (21-22), (24-25), (26-27-28); 

E: Te (4,11,12,13). 



* According to the original results of Kriiger et al. (1994) (i.e. /v = 1; see Sect. 2.2.1); in all other models /v = 0.588 is adopted. 



to H-ionising photons (/i(H+) < 1), and presents line ra- 
tios [Oiii] A5007/HeiiA4686 and He ii A4686/H/3 that at 
late times are very different from those predicted for the 
PN model of Fig. These differences are mainly due 
to the fact that, in comparison to the H-burning central 
star, the maximum flux of ionising photons emitted by the 
He-burning star is attained later (because of the longer 
evolutionary time-scales), when the nebular material has 
expanded to larger radii. 

3.2. Detection probabilities of different PNe 

The relative probability of observing PNe of different pro- 
genitors (i.e. stars of different masses, ages, and metallic- 
ities) constitutes a basic aspect in the interpretation of 
PNe data. As briefly mentioned in Sect. H, the highest 
probabilities of observing the PN stage should be reached 



when we have a suitable combination between the evo- 
lutionary time-scale of the central star, and the dynam- 
ical time-scales of the expanding nebulae: A fine tuning 
between the two can increase either the lifetime or the 
surface brightness of a PN, both factors enhancing the 
detection probability. 

Actually, the problem of determining a "detection 
probability" is much more complex. For instance, in the 
case of distant (point-like) nebulae, the detection proba- 
bility should correlate with the absolute luminosity, rather 
than with the surface brightness. Also the choice of the de- 
tection method is crucial in determining which PNe will 
be better sampled by a given survey. Since our models pre- 
dict a lot of different properties of PNe, we are also able 
to face this question in a more generic way: we can con- 
struct synthetic samples of PNe containing all their basic 
observable properties, which are then be filtered adopting 
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Fig. 11. Predicted evolution of basic parameters corresponding to a PN model ionised by a H-burner central star with 
mass of ~ 0.6 Mq (model 4). See Sect. ^ for a description 



some given observational criteria. Developing such a tool, 
however, is beyond the scope of this paper. 

Anyway, it is worth getting a first insight into the prob- 
lem. To this aim, in Fig. |l^ we plot the evolution of two 
PNe models in the HR diagram, corresponding to H- and 
He-burning tracks (models 4 and 9, respectively). Along 
the tracks, we mark the CSPN position at equally spaced 
intervals of age (700 yr) by means of circles whose area is 
proportional either to the emitted line luminosity L\ (left 
panels) or to the PN average surface brightness S\ (right 
panels) in a given emission line (either H/3 or [O ill] A5007). 
S\ is simply assumed to be proportional to the quantity 
Lx/Rf^j^, where i?ion is the total nebular radius. 

The figure gives a clear indication of the regions of 
the HR diagram where these PNe are more likely to be 
found, as well as of their relative luminosities and sur- 
face brightnesses. In the case of the H-burning track, the 
probability of detection in the H/3 line should be more or 
less evenly distributed in effective temperature, as long 
as log Toff 4.2 and the CSPN is crossing the HR dia- 
gram horizontally: In this part of its evolution, the nebula 



evolves at slightly increasing Lh,3, and slightly decreasing 
S'h/3. Also, there is a non- negligible probability of finding 
these PNe at later ages, when their central stars reach 
the initial part of the white dwarf cooling sequence (at 
logL/L© ~ 2.3, log Toff ~ 5.05) and the evolutionary 
speed slow down. However, in this part of the HR dia- 
gram, S'h/3 is extremely low. 

The same H-burning nebula should be much brighter 
in the [O ill] A5007 line. Remarkably, during the horizontal 
crossing of the HR diagram, the detection of [O ill] would 
be favoured in a more limited temperature interval, i.e. 
5.0 ^ logTeff ^ 4.2; at higher temperatures, the [O ill] 
line temporarily disappears. This effect is caused by the 
photoionisation of 0++ into and is explained in 

more detail in Sect. 4.(;. 



Figure also illustrates the striking differences be- 
tween the H- and He-burning cases. The PN with the He- 
burning track have line luminosities that are just slightly 
lower than its H-burning counterpart. However, except 
for the very initial stages, the He-burning model presents 
much lower surface brightnesses. This is caused by the 
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Fig. 12. The same as in Fig. but for a He-burner central star with mass of 0.6 (model 9) 



longer evolutionary time-scales of the He-burning central 
stars. Moreover, the He-burning CSPN makes a long pause 
in its evolution amid the horizontal crossing of the HR di- 
agram, at Teff ~ 4.7, period in which the nebula should 
have a very high probability of being detected. During this 
pause, the line luminosity is almost constant, but the sur- 
face brightness decreases a lot, due to nebular expansion. 
Thereafter, the CSPN resumes its evolution at a slower 
pace than in the H-burning case. In summary, PN with 
He-burning tracks should present a very odd distribution 
in the HR diagram: most of them should concentrate in the 
logTcff ^ 4.7 region, with a small fraction of high-surface 
brightness objects appearing at 4.7 ^ logT^ff ^ 4.3. This 
trend should be present in both H/3 and [O ill] A5007 hues. 
However, similarly to the H-burning case, the [O ill] line 
disappears in the upper-left section of the HR diagram, 
when Tcff Z 4.9 and logL/Lg ^ 3. 



4. Comparison with observations 

The next sections are dedicated to present the results of 
PNc model calculations in comparison to the observed 
data, which are used to both calibrate the free param- 
eters and, more generally, test the theoretical prescrip- 
tions. As to the calibration of the parameters, we restrict 
the comparison here to a sample of Galactic PNe, so as 
to minimise possible population effects (due e.g. to dif- 
ferent mean metallicity) that may be present among PNe 
samples of nearby galaxies. 



4.1. The ionised mass - nebular radius relationship 

PNe display a positive correlation between the ionised 
mass and the nebular size, which is shown in Fig. ^ for 
Galactic PNe. Such relation offers also a method to esti- 
mate the distance of PNe, once the ionised mass and the 
angular radius can be measured (Maciel & Pottasch 1980; 
see also Kwok 1985). The mass-radius relationship is nat- 
urally explained in the framework of the interacting-winds 



p. Marigo et al.: Evolution of Planetary Nebulae I. An improved synthetic model 



19 




log Ttrf, log T^^ff 

Fig. 13. Representation of the relative probabilities to detect a PN observed either in H/3 (top panel) or in the [O 
III] A5007 line (bottom panel), as its central star evolves along its post-AGB track, assumed either H-burning (model 
4) or He-burning (model 9, that corresponds to the least luminous track in each panel). Circles are drawn at equally- 
spaced time intervals of 700 yr, their area being proportional to the current values of the luminosity (left panels), and 
surface brightness (right panels) in each emission line 
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In this study we use the observed relation first to derive 
some calibration of ^a^Ib' ^'^^ which we presented three 
possible prescriptions in Sect. 2.2.1. To this aim, we cal- 



0.01 0.1 1 

Fig. 14. Ionised mass - nebular radius relationship. 
Observed data for Galactic PNe are taken from the compi- 
lation by Zhang (1995; squares), and the sample presented 
by Boffi. & Stanghellini (1994; triangles). Model predic- 
tions are displayed for various prescriptions for and 
itr (see Table ^) 



model (VK85), which couples the dynamical expansion of 
the nebula with the ionising properties of the central star. 



culate PN models all corresponding to a CSPN with mass 
of about 0.6 Mq, since this should be representative of the 
typical PN population in the Galactic disk (see Stasihska 
et al. 1997, and references therein). The resulting values 
for l-AGB i) const. = 15 km s~^ (models 1 and 5), ii) 
F{Z) ^ 27 km s~^ (model 2) according to Bressan et al. 
(1998), with Z being the effective metal abundance (higher 
than the initial value because of dredge- up episodes); iii) 
F{L) ^ 18 km s^^ (model 3) according to Kriiger et al. 
(1994V We can notice that all the corresponding curves in 
Fig. |lj lie below the mean line of the observed data, even 
when a shorter transition time (ttr = 500 yr, model 5) is 
adopted to reduce the discrepancy. 

We find that the mean mass-radius relation can be 
well reproduced by decreasing V^^, e.g. applying a mul- 
tiplicative factor /v — 0.588 to the original Kriiger et al. 
data (so that Vagb ~ 10 model 4 ; see also 

Fig. p. Most of our PN model calculations, presented in 
the next sections, are carried out by adopting this cali- 
brated prescription. An additional advantage, compared 
to the case of a constant velocity, is that such a choice 
produces a natural dispersion of the expansion velocities 
at the end of AGB, which seems more consistent with the 
observed scatter of the expansion velocities of young PNe 
at small radii (see Fig. ^l|) . 

With this calibration for , we show in Fig. ^5] the 

predictions for several choices of the other model param- 
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eters (quoted in Table g). In general, all tracks exhibit 
a common feature in the Mion--Rion plane, i.e. a change 
in the slope after an initial steep increase of the ionised 
mass. This point marks the transition from ionisation- to 
density-bounded shell configurations. In other words, dur- 
ing the early epochs the ionisation front proceeds forward 
but keeping behind the outer shock, so that the nebula is 
optically thick to the ionising photons, and the outermost 
part of the shocked shell is still neutral. 

As soon as the ionisation front reaches the outer rim of 
the shell, the unperturbed (low-density) AGB wind is also 
ionised and the nebula becomes optically thin to the con- 
tinuum photons emitted by the central star. During these 
stages the whole shell is ionised, and the theoretical mass- 
radius tracks merely represent the evolution of the shell 
parameters. The almost linear behaviour proceeds until 
when - due to the decrease of the ionising flux from the 
central star - the ionised mass starts to decrease as well, 
and correspondingly the ionisation front recedes inside the 
shell. At later ages, Mjon increases again as the nebula ex- 
pands and its density decreases. Finally, when the AGB 
wind is entirely engulfed by the shocked shell, Mjon re- 
mains constant at increasing radius (horizontal part of 
the tracks). 

The evolution of the optical properties of PNe along a 
mass-radius track is also shown in Fig. ^ together with 
the observed data for the sample of Galactic PNe studied 



by Mendez et al. (1992), who estimate the /i-factor char- 
acterising each object. Both the location and the slope of 
the theoretical relationship for optically thick/thin sources 
are in excellent agreement with the observed data. 

Clearly, this general trend may be affected by vari- 
ous parameters. Let us first discuss the models illustrated 
in panel A of Fig. |l5|, all corresponding to a central star 
with mass of « O.QMq. Most of these PNe models ac- 
count for the mean observed mass-radius relation quite 
well, which is an expected result as they should be rep- 
resentative of the typical Galactic PNe. For a given tran- 
sition time (itr = 1500 yr), the largest differences show 
up at the lowest and largest values of Mjon, whereas the 
models almost converge during the optically-thin stages. 
For a significantly longer ttr (8000 yr) the observed data 
are intercepted by the theoretical track (model 8) only at 
the largest i?ion, since the nebula starts to be ionised when 
it has already expanded to relatively large dimensions. 

The dependence on the mass of the central star can 
be visualised in panel B of Fig. 15. It turns out that the 



radius at which ionisation first sets in does not vary much 
with McsPN, being essentially determined by the transi- 
tion time (the same in all models) . A clear effect of Mcspn 
is instead related to the part of the tracks following the 
transition from optically-thin to optically-thick configu- 
rations (the first relative maximum in the curves). At in- 
creasing stellar mass, this transition occurs at shorter radii 
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Fig. 16. Ionised mass as a function of the nebular radius 
for the sample of Galactic PNe analysed by Mendez et 
al. (1992). Optically-thin/thick sources are distinguished 
by means of empty/filled symbols, respectively. Along the 
theoretical track, corresponding to model 18, the expected 
optical properties are indicated (thick line: optically-thick 
configuration with /i = 1; thin line: optically-thin config- 
urations with < 1 ) 
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Fig. 17. Relative shell thickness as a function of the 
nebular radius. Measured radii are taken from Boffi & 
StanghcUini (1994; triangles); Zhang (1995; squares); and 
Gorny et al. (1997; circles). Corresponding semi-empirical 
determinations of 6R are those derived by Zhang & Kwok 
(1998). A few theoretical predictions are shown for com- 
parison (see Table H) 



reflecting the faster evolutionary speed. Moreover, we can 
note that the non-monotonic behaviour of the tracks can 
explain part of the observed scatter of Mjon at given i?ion ■ 

The effect of the transition time can be analysed in 
panel C. In this panel, we can notice that the observed 
PNe with the lowest masses and radii would be consistent 
with relatively massive central stars and quite short tran- 
sition times {ttr = 500 yr; models 22, 29). Models with 
lower McsPN and much longer ttr (10000 yr; models 15, 
14) succeed in explaining the few data at quite large -Rion 
and relatively low Afjon- 

Finally, we plot in panel D the expected mass-radius 
tracks for He-burning post-AGB tracks. These models are 
consistent with the observed data as well as those with H- 
burning central stars. Comparing models 21 (H-burner) 
and 23 (He-burner) - which are calculated with the same 
set of the other parameters - we can see that, contrary to 
the H-burning case, in the He-burning model the ionised 
mass always increases with radius and the nebular shell, 
once becomes optically thin, remains density-bounded up 
to the largest radii. This can be explained considering that 
He-burning tracks are, on average, characterised by quite 
longer evolutionary time-scales, so that the drop of the 
ionising flux emitted by the central star occurs when the 
expanding shell has already reached large radii. 



4.2. The nebular shell thickness 

PNe exhibit a large variety of morphological and geometri- 
cal features (e.g. Schwarz et al. 1992), that clearly cannot 
be handled by models that assume spherical symmetry. 
Nevertheless, we can at least test our predictions concern- 
ing an important geometrical parameter, namely the rel- 
ative shell thickness SR (deflned by Eq. ^). As already 
described in Sect. 2.7, our model considers the effect of 



ionisation in determining the thickness of the shell. 

To construct Fig. |l^, we combine determinations of 
nebular radii taken from different sources with the re- 
cent semi-empirical determinations of SR (for the primary 
shell) obtained by Zhang & Kwok (1998) in their mor- 
phological analysis of a sample of Galactic PNe. For each 
of them, Zhang & Kwok (1998) derive relevant morpho- 
logical parameters, by constraining simulated optical and 
radio images to reproduce the observed ones. It turns out 
that SR mostly varies over the range [0.1, 0.8], with a dis- 
tribution peaked at around 0.45. 

We can see that the theoretical predictions (a discus- 
sion of the trends is given in Sect. 2.7) account well for 
the location of the observed data in the SR-Rion plane. 
Together with the good reproduction of the ionised mass- 
radius relationship, this fact implies that both the dynam- 
ical and ionisation properties of PNe should be reasonably 
well represented in our model. 
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Fig. 18. Electron densities as a function of the nebu- 
lar size. Observed data for Galactic PNe are taken from 
the compilation by Phillips (1998). Model predictions are 
shown for comparison. The straight line, Nc oc R^"^, is 
plotted to indicate the expected slope for constant nebu- 
lar mass 
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Fig. 19. Zanstra discrepancy between Tz(H) and 
rz(Heii). Observed data for Galactic PNe are taken from 
de Freitas Pacheco et al. (1986; triangles); Gleizes et 
al. (1989; starred symbols); and Mendez et al. (1992; 
squares). Model predictions for different choices of the pa- 
rameters are shown for comparison 



4.3. Electron densities 

Figure |l^ displays the radial variation of the nebular elec- 
tron densities. Empirical determinations of Nc are all de- 
rived from the [S ii] A6717/6730 forbidden line ratio, which 
should sample the outer nebular regions. In this way, we 
expect that the considered PNe sample is little subject to 
the effects of radial stratification in ionisation and den- 
sity - present when different lines are used - as well as 
not dependent on the nebula filling factor, which is the 
case when is estimated from the radio continuum/H/3 
fluxes. 

From the observed distribution of the data points it 
should be noticed that, first of all, the relation between 
and i?ion shows a relatively large dispersion, which seems 
to reduce towards larger radii. Second, as pointed out by 
Phillips (1998), there is clear evidence for a change in the 
mean slope of the radial gradient at i?ion « 0.1 pc. Both 
features are quite well reproduced by our models, showing 
the evolution of the average density of the ionised shell as 
function of the radial size of its outer boundary. The data 
points located at i?ion < 0.1 pc would be consistent with 
a mixture of both optically-thick and thin nebulae, the 
latter being characterised by a more gradual decrease of 
Ne with radius. Indeed, this feature would explain the ob- 
served scatter of the points at smaller i?ion- For Rio^ > 0-1, 
the reduced dispersion of the observed data which, on av- 
erage, seem to follow the relation Nc cx could be 
interpreted by the gradual convergence of all PNe to the 



condition Mion = Msh = AMagb, i-e. the shocked ionised 
shell embraces all the AGB ejecta, so that its mass remains 
constant. Finally, we note that the theoretical Nc — i?ion 
tracks show a systematic trend with the mass of the cen- 
tral star, i.e. at smaller radii the largest densities are at- 
tained for more massive CSPN. However, it should be re- 
marked that the predicted N^ refer to average electron 
densities, which do not account for intrinsic density gra- 
dients inside the nebulae. 

4.4. PN central stars: temperatures and luminosities 

Determining the main parameters of the central stars of 
PNe , i.e. effective temperature and luminosity, means to 
cast light on important aspects of both nebular and stellar 
evolution. For instance, by assigning a position on the H- 
R diagram to the ionising star, it is possible to derive 
the stellar mass - from the comparison with theoretical 
post-AGB tracks - and estimate (an lower limit to) the 
PN dynamical age, if the expansion velocity and nebular 
sizes are measured. The combined information on stellar 
mass and PN age, offers also the possibility to discriminate 
between H- and He-burning stars (e.g. Dopita et al. 1996). 

Unfortunately, the determination of luminosities of 
Galactic CSPN is indeed problematic, mostly due to 
the large uncertainties that often affect the observational 
methods to estimate effective temperatures and distances. 
The effective temperature of a CSPN can be determined 
either directly, by fitting observed spectra with model at- 
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Fig. 20. Zanstra temperatures 
(Tz(H) in the left panels and 
Tz (He ii) in the right ones) and 
luminosities (/ini) of PN central 
stars. The upper panels present 
the nebulae with CSPN masses 
0.65 < McspN < 0.7, whereas the 
lower panels refer to those with 
0.55 < McsPN < 0.6. Observed 
data points are taken from Mendez 
et al. (1992), and shown with 
filled squares for > 0.7 and 
empty squares for /ijj < 0.7. 
PNe of low-excitation (with Tz(H 
only) and high-excitation (with 
both rz(H and Tz(Heii)) are 
denoted with squares and circles, 
respectively. Model predictions 
are superimposed distinguishing 
optically-thin / thick configurations 
(thin/thick solid lines) in the H 
and He+ continua, together with 
the corresponding logL — logTcff 
tracks (dotted line). The masses 
of the central stars are chosen 
consistently with those estimated 
by Mendez et al. (1992) for the 
sample PNe 



mospheres, or with the aid of indirect methods, such as 
the Zanstra and energy-balance method (see, e.g., Gathier 
& Pottasch 1989). The former, first introduced by Zanstra 
(1931), is based on the the assumption that the nebula is 
optically thick to the ionising continuum under consider- 
ation - the H and He ii Zanstra temperatures Tz (H) and 
Tz (He ii) are usually determined - and requires the mea- 
surement of both nebular and central star optical emis- 
sions. The latter, first discussed by Stoy (1933), is based 
on the assumption that the ratio of the total energy emit- 
ted by a nebula in forbidden lines to the energy emitted in 
one of the hydrogen recombination lines depends on the 
temperature of the central star. Then, once Teff is esti- 
mated, the luminosity of the CSPN can be derived, for 
instance, i) from the apparent magnitude at a given wave- 
length provided that the distance is known, ii) from the 
comparison with post-AGB tracks in the log 5 — log Teff di- 
agram if estimates of the surface gravity g and mass of the 
CSPN are available; iii) by equating the nebular emission 
over the whole spectral range to the radiation emission 
from the central star, with the assumption that the neb- 
ula is ionisation-bounded (see, e.g. Gathier & Pottasch 
1989, Mendez et al. 1992 for description and application 
of the methods). 

The so-called "Zanstra discrepancy" between Tz(H) 
and Tz(Heii) - i.e. usually Tz(H) < Tz(Heii) - has been 
a matter of debate for long time (see, e.g., Stasihska & 



Tylenda 1986; Gathier & Pottasch 1989; Mendez et al. 
1992; Gruenwald & Viegas 2000 and references therein). 
The discussion aims at identifying the main causes of such 
discrepancy and their role, namely: i) the condition of op- 
tical thinness to the H continuum; ii) the possible excess 
of photons able to ionise He+ - with energy hv > 54.4 eV 
- emitted by the CSPN, compared to the black-body spec- 
tral energy distribution; iii) the differential dust absorp- 
tion in the nebula. 

Figure |l^ illustrates the Zanstra discrepancy for 
Galactic PNe. Model predictions for rz(H) and Tz(Heii) 
are derived from the conditions: 

A.hQ(H°) = iTT r ' :?4^d(/.^) (44) 

MHe.Q(He+) ^ 4. r ^-[^f^^")] d(H (45) 

where is the Planck function and the other quanti- 
ties have the usual meaning. It should be noted that we 
do not strictly follow the same procedure adopted by ob- 
servers to obtain the Zanstra temperatures, which usually 
involves the ratio between nebular and stellar emission at 
particular wavelengths (the so-called Zanstra ratios; see 
e.g. Mendez et al. 1992). However, Eqs. (0) and (^) share 
the same basic assumptions, i.e. the nebula is considered 
in equilibrium state between ionisation and recombination 
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processes, and optically thick to the ionising continuum 
(ionisation-bounded) . 

The theoretical tracks in Fig. |l^ account for the ob- 
served trend of the Zanstra discrepancy. From the analysis 
of the models it is worth reporting that i) at lower Tz(H) 
the Zanstra discrepancy increases with both the stellar 
mass and the transition time, and ii) the highest values 
of Tz(Heii) are attained by the most massive models; iii) 
there is some hint for an inverse Zanstra discrepancy (i.e. 
Tz(H)> Tz(He ii)) which may occasionally occur when the 
central star overcomes the hottest knee of its post-AGB 
track and starts to decline in luminosity. 

Figure ^ contributes to investigate the causes of the 
Zanstra discrepancy. To this aim, we adopt the sample 
of 23 Galactic PNe analysed by Mendez et al. (1992), for 
which a large number of parameters are determined, in- 
cluding: the nebular masses and radii (see also Sect. |4.lD , 
the masses of the CSPN, the effective temperatures, the 
Zanstra temperatures, etc. The observed data for T'z(H) 
and Tz(Heii) are plotted as a function of the quantity 
yLa as defined by Mendez et al. (1992; see their table 4), 
which is equal to /xni according to their equation 12. The 
comparison with the theoretical predictions is then made 
for representative models whose stellar masses fall within 
the range indicated by Mendez et al. 

The general accordance between the observed data and 
their expected location is remarkably good. The observed 
dispersion of the data in the log(/iH^) — logTz(H) plane 
(see top-left panel of Fig. |2^) is naturally accounted for by 
the models, which exhibit notable deviations from the hor- 
izontal parts (at constant luminosity) of the log L — log Tog 
post-AGB tracks. Such deviations show up as soon as the 
ionised shell becomes optically thin (i.e. /j,h < 1) to the H 
continuum, which confirms the conclusions already drawn 
by Mendez et al. (1992). Note, however, that there is some 
offset in Teft between observations and models, i.e. the lat- 
ters get optically thin at lower effective temperatures than 
indicated by the data. 

In the log(/ZH-^) — log 7z (Hen) plane (see top-right 
panel of Fig. ^0|) the models succeed in reproducing the 
corresponding data which, in this case, are expected to 
correspond to optically thick configurations to the He'*" 
continuum. It is also worth noticing that the agreement 
holds as well when the data are analysed as a function of 
the mass of the CSPN, as is clear by comparing the top 
panels (hig her masses) to the bottom ones (lower masses) 
of Fig. |20|. Finally, a general convergence is found with 
respect to the discrimination between optically thin/thick 
nebulae and their location in the observational planes. 

4.5. Expansion velocities 

The expansion velocities of PNe - usually derived from 
emission lines like Hi, Hell, [Sii], [On], [Ariv], [Neiii], 
and more often [O ill] - are found to range from few to 
40 - 45 km s~^, without any well-defined correlation with 
the nebular size (see Fig. 21). 



Moreover, velocity estimates based on lines produced 
by ions with different ionisation potentials, have revealed 
the existence of velocity gradients across the nebulae 
(Gesicki et al. 1998). In general, the velocity is found to 
increase with the ionisation potential of the emitting ion, 
i.e. going from the inner to the outermost layers of the 
nebula. The fact that a large scatter (shown in Fig. |l]) 
still characterises the velocity data, though derived from 
the same lines, likely points to a relatively broad range of 
parameters affecting the dynamics of PNe. 

We attempt to investigate this point on the basis of the 
analytical prescriptions already presented in Sects. 2.4 and 
p.7| . We just recall that, although the adopted dynamical 
description is a simplified approach that cannot account 
for the complex velocity structure shown by hydrody- 
namical simulations (e.g. Schmidt- Voigt & Koppen 1987b; 
Marten & Schonberner 1991), nevertheless it should pro- 
vide a reasonable representation of the average dynami- 
cal properties of the expanding nebulae. Furthermore, a 
clear improvement over similar analytical models is that 
the expansion of the nebula is not expressed with a sin- 
gle velocity, but the velocities of both the inner and outer 
rims of the ionised shell (Vg*, and Vion respectively) are 
predicted. In this way, the analytical model keeps track, 
to some extent, of the real velocity gradient inside the 
nebula. 

An example is shown in Fig. |^ (top panel), where 
the two rim velocities are plotted as a function of -Rion- 
We can immediately notice the relatively large velocity 
range comprised by the two curves. In particular, quite 
low velocities (few km s^^) are attained by the contact 
discontinuity during the early post-AGB stages, when the 
increased pressure due by ionisation slows down its ex- 
pansion (see Sect. 2/7 for more discussion). Conversely, 
the ionisation front is always accelerating outward. 

On the other hand, for purpose of easier comparison 
with observed data, we may also define an average expan- 
sion velocity 



(14xp) = 



47r(i?i: 



ion 



Rf) 



V{r) Airr'^dr 



(46) 



following Schmidt- Voigt & Koppen (1987b). In practice, 
it is a weighted mean over the nebular (emitting) volume, 
under the assumption that the velocity varies linearly be- 
tween the two extrema Vg and VJon- According to this 
definition (I4xp) is closer to VJon, that is the maximum 
velocity. 

The dependence of (Vcxp) on the GSPN mass (for H- 
burning models) is evident already during the initial stages 
of the nebular expansion, so that the highest velocities are 
reached by the most massive models (see middle panel of 
Fig. ^l|). We can also notice that all the rising branches 
of the curves concentrate on a narrow range in radius, 
just refiecting the fact that the transition times of these 
models are comparable (i.e. 1500 yr in most cases). After 
the initial rise and peak, the curves start then to decline 
and finally flatten out at velocities of about 20 - 30 km 
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Fig. 21. Expansion velocities as 
a function of the nebular size of 
Galactic PNe. Data (circles) de- 
rived from the forbidden lines of 
[O III] and [O ii] are taken from the 
compilation by Weinberger (1989). 
Model predictions are shown for 
comparison for several choices of 
the parameters. See text for further 
details 



s^^, which coincides with the mean range of the observed 
data. However, we remark again that this is a just a mean 
trend. For the same set of models the velocities of the 
outer rims are shown for comparison (dotted lines). 

The initial radius of each velocity curve essentially rep- 
resents the stage at which the nebula starts being ionised 
by the central star. This radius, i?ion,Oi can be expressed in 
first approximation as 1^X36 ^ ^tr, i-e. the product of the 
velocity of the ejected AGB material and the transition 
time. At increasing transition times, ionisation sets in at 
larger radii (compare model 4 with 8, having tn = 1500 
and 8000 yr, respectively). The same effect is produced 
for lower values of the ejection velocity, as we can see in 
the bottom panel by comparing model 14 (solid line) with 
model 15 (short-dashed line), having both tn = 10 000 yr, 
but VaxI'b ~ 6, 17 km s~^ respectively. It follows that at 
least a part of the observed velocity spread of PNe at any 
radius can reflect differences in the evolutionary speed of 
the central star, as well as be the result of differences in 
the ejection velocity on the AGB. 

It is worth also noticing the different behaviour of 
(Vexp) of PNe ionised by central stars that are either H- 
(e.g. model 4) or He-burners (e.g. model 9). On average, 
we expect that at given nebular size a PN ionised by a 
He-burning CS has a lower expansion velocity than that 
surrounding a H-burning CS of the same mass. This in 
turn would translate into longer dynamical ages, - de- 



fined as tdyn — ^ion/Kjxp " of PNc with He-burning CSs. 
This age difference offers a viable way to discriminate be- 
tween the two classes of stars (see Dopita et al. 1996 for 
an application to PNe in the LMC). 

On the basis of all mentioned points it is clear that, 
even if we are interested in average dynamical properties 
of PNe, the often adopted assumption of a constant ex- 
pansion velocity can lead to substantial errors in the es- 
timation of the true dynamical ages. Actually, it will be 
interesting to use our model results in a future analysis 
to address the discrepancy between the evolutionary and 
dynamical ages of PNe (see e.g. McCarthy et al. 1990). 

4.6. Nebular line ratios 

The ionisation structure and line emissivities are among 
the most basic properties of our PNe models, and can be 
easily checked by comparison with observations. 

4.6.1. The Heii/H/3 vs. [Oiii]/Heii anticorrelation 

Let us start with Figure |2^, which combines the in- 
tensities of the three optical emission lines here consid- 
ered, i.e. recombination (H/3, HeiiA4686) and forbidden 
([Oiii] A5007) lines. 
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Fig. 22. Nebular line ratios in Galactic PNe. Data are taken from the Strasbourg-ESO Catalogue of Galactic Planetary 
Nebulae (Acker et al. 1992). Theoretical tracks are superimposed for comparison. Let panel: models 4-8-9-14-28, all 
calculated assuming a constant electron temperature To = 10 000 K. Right panel: models 4-11-12-13, calculated with 
Tc — 8 000, 10 000, 12 000 and 13 000 K (from higher to lower tracks). The arrow indicates the direction of increasing 
Tcff 



First of all, we notice that the observed data do not 
present a large scatter, but rather follow a well-defined 
relation. There is an evident anticorrelation between the 
He II A4686 /H/3 and [O iii] A5007/He ii A4686 intensity hue 
ratios, until the former becomes almost constant 1) for 
low values of [O iii] A5007/Heii A4686. 

The interpretation of such trend can be easily de- 
rived when considering that the ionisation potentials of 
He+ (54.4 eV) and 0++ (54.9 eV) are almost identi- 
cal. It follows that the Stromgren volumes F(He+) and 
V{0~^'^) nearly coincides and the lines HeiiA4686 and 
[O III] A5007 are mainly produced within separate nebu- 
lar regions, whose complementarity in emission volumes 
(i.e. when the innermost ^(He"'""'') increases, the outer 
V(0"'""'") decreases and vice- versa) explains the observed 
anticorrelation between the respective line intensities. 

It is also clear that since 116+"^ is present in the neb- 
ula only if the central star is sufficiently hot, the observed 
trend in Fig. ^ implicitly contains a temperature scale 
(increasing along the direction of the arrow) , and as a mat- 
ter of fact, it has been used as a semi-empirical method to 
determine the temperature of high-excitation CSPN (see 
e.g. Gurzadyan 1988). 

It is really striking how the empirical relation is so 
well reproduced by our models, which stands on few sim- 
ple prescriptions for the ionisation and emission nebular 
structure (see Sect. 2.9). In this respect we draw the at- 
tention to the following points: 

The flattening of the tracks towards the lowest 
[O III] A5007/Heii A4686 values occurs when the nebulae 



become optically-thin to both the H and IIe+ continua. In 
fact, as the emission volumes for these ions are the same, 
the He ii A4686/H/3 intensity ratio depends only on atomic 
parameters (i.e. recombination coefficients and chemical 
abundances). 

it is ev- 



Finally, comparing the two panels of Fig. 22 



ident that for given He ii A4686/H/3 intensity ratio, the 
relative intensity of the forbidden line is a strong func- 
tion of the electron temperature (right-panel), whereas it 
is almost independent of other parameters (e.g. transition 
time, mass of the CSPN, H-/He-burning conditions, etc.). 
As expected, most of the observed data can be reproduced 
if Tc is varied from about 8 000 to 13 000 K, in agreement 
with the typical temperature range indicated by observa- 
tions. 



4.6.2. The evolution of the /(He ii A4686)//(H/3) line 
intensity ratio 

Figure H shows the behaviour of the /(He ii A4686) //(H/3) 
intensity ratio as a function of the nebular radius. At 
smaller radii (i?ion ^0.1 pc), the observed data scatter 
over a large range, levelling up to values of 1. This trend 
is explained by the models as the result of the progres- 
sive increase of the stellar temperature at nearly constant 
hmrinosity, so that for Tos > 60 000 K helium starts to 
be doubly ionised and the recombination line He ii A4686 
starts to be detectable. 

At this stage the nebula is already optically thin to 
the H-contimmm (see Fig. |l^). As the temperature in- 
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Fig. 23. Intensity line ratio as a function of the nebular 
size. Observed data refer to the sample of Galactic PNe 
analysed by Gorny et al. (1997). Model predictions are 
shown for different choices of the CSPN mass 

creases, the Stromgren volume ^(He "'""'") grows (and the 
HellA4686 line intensity strengthens) until when it em- 
braces the whole nebula, which then becomes optically 
thin also to the He-continuum. This occurs in correspon- 
dence to the maximum temperature attained close to the 
knee of the post-AGB track. As already mentioned, this 
condition yields the /(He ii A4686)//(H/3) line ratio to be 
practically constant and close to unity (as shown by the 
horizontal part of the tracks). 

As soon as the stellar temperature starts declining, 
the Stromgren radius i?ion(Heii) recedes inward, the 
He II A4686 intensity weakens the H/3 line is still relatively 
strong (and the nebula still optically thin). We can no- 
tice that at some point in all models the decrease of the 
/(He II A4686) //(H/?) line ratio slows down while the neb- 
ular radius is still increasing. This happens when the H/3 
line intensity is getting weaker as well and the nebula re- 
turns to be optically thick to the H-contimium. Such pre- 
diction is in perfect agreement with the observed lack of 
very low /(Hell A4686)//(H/3) line ratios at larger radii 
(i.e. Rion > 0.2 pc). 

5. Concluding remarks 

In the present work we have presented in detail a new 
synthetic model for FN evolution. On the one hand, this 
model is kept simple by using convenient approximations 
for most physical processes, and by limiting their descrip- 
tion to analytical relations that can easily be integrated 
with time. On the other hand, the model improves upon 
previous synthetic PNe models in the attempt i) to consis- 



tently couple the FN phase to the previous AGB evolution 
of the stellar precursor, and ii) to account for dynamical 
effects (like that produced by ionisation) that are believed 
important to the evolution of PNe. 

Although the present description considers many as- 
pects of the FN evolution, it cannot, of course, be consid- 
ered exhaustive. Some approximations - like the assump- 
tions spherical symmetry and constant electron tempera- 
ture, and the somewhat simple description of the dynam- 
ical wind interactions - are in fact too crude if compared 
with the results of detailed hydro dynamical and photoion- 
isation codes. However, this is the due price for fast com- 
putations. We remark that such synthetic codes proba- 
bly provide the only way, at present times, to avoid the 
remarkable time requirements of hydro dynamical simula- 
tions. 

This compromise is similar to that in the field of AGB 
evolution (see e.g. Groenewegen & Marigo 2001): AGB 
stars can be described either by i) complete stellar evo- 
lution codes, that due to their demand of CPU-time can 
be run only for few choices of stellar parameters, or by ii) 
synthetic AGB codes, that treat ill-unknown quantities by 
means of a parametric approach, and that can be easily 
run for any choice of parameters. Synthetic AGB mod- 
els, although always recognised as crude tools, have had a 
crucial role in the understanding of many aspects of AGB 
evolution, and are widely considered to be fundamental in 
the field. 

Similarly, our synthetic PNe code is fast enough to al- 
low the systematic computation of extended grids of PNe 
tracks - for instance, for different choices of mass and 
metallicity distributions. The few free parameters (tran- 
sition time, electron temperature, H-/He-burning tracks, 
etc.) can be constrained within reasonable limits by com- 
paring the results with observations, in a way similar to 
the one we have presented in this paper. Anyway, none of 
these parameters is expected to be fixed for a given PNe 
population. In fact, the most promising way to explain 
the wealth of PNe characteristics, is that of performing 
Monte-Carlo simulations in which the free parameters are 
given reasonable probability distributions. 

In the forthcoming papers, we will apply this approach 
to produce synthetic populations of PNe in different galax- 
ies. Main targets will be: 

— To further constrain the models by means of more 
detailed comparison with the observed properties of 
Galactic, Magellanic Cloud, and other PNe popula- 
tions. 

— To explore the expected distribution of chemical abun- 
dances among PNe, that may help to constrain syn- 
thetic TP- AGB models and their chemical yields. 

— To simulate the PNe luminosity function in different 
galaxies, in order to test their suitability as extragalac- 
tic distance indicators. 
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